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We study two weighted graph coloring problems, in which one assigns q colors 
to the vertices of a graph such that adjacent vertices have different colors, with a 
vertex weighting w that either disfavors or favors a given color. We exhibit a weighted 
chromatic polynomial Ph(G, q, w) associated with this problem that generalizes the 
chromatic polynomial P(G,q). General properties of this polynomial are proved, 
and illustrative calculations for various families of graphs are presented. We show 
that the weighted chromatic polynomial is able to distinguish between certain graphs 
that yield the same chromatic polynomial. We give a general structural formula for 
Ph(G, q, w) for lattice strip graphs G with periodic longitudinal boundary conditions. 
The zeros of Ph(G, q, w) in the q and w planes and their accumulation sets in the 
limit of infinitely many vertices of G are analyzed. Finally, some related weighted 
graph coloring problems are mentioned. 



I. INTRODUCTION 

Recently we have formulated two weighted graph coloring problems in which one assigns q 
colors to the vertices of a graph such that adjacent vertices (i.e., vertices connected by an edge 
of the graph) have different colors, with a vertex weighting w that either disfavors (for < 
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w < 1) or favors (for w > 1) a given color We label these with the abbreviations DFCP 
and FCP for disfavored-color and favored-color weighted graph vertex coloring problems. 
Since all of the colors are, a priori, equivalent, it does not matter which color one takes to 
be subject to the weighting. In the present paper we study these interesting weighted graph 
coloring problems in detail. An assignment of q colors to the vertices of a graph G, such 
that adjacent vertices have different colors, is called a "proper g-coloring" of the vertices 
of G. We analyze the properties of an associated weighted chromatic polynomial that we 
denote Ph(G,q,w), which generalizes the chromatic polynomial P(G,q) and constitutes a 
w-dependent measure, extended from the integers to the real numbers, of the number of 
proper g-colorings of the vertices of G. In the weighted graph coloring problem, with q G N + 
being the number of colors, for a given graph G, P(G, q) is a map from N + to N, while 
Ph(G,q,w) is a map from N + x J to 1, where I denotes the DFCP interval < w < 1 
or the FCP interval w > 1. In both cases, one can formally extend the domain of q and 
w to R or, indeed, C, and the latter extension is, in fact, necessary when one analyzes the 
zeros of P(G,q) or Ph(G,q,w), respectively. The polynomial Ph(G,q,w) is equivalent to 
the partition function of the g-state Potts antiferromagnet on the graph G in an external 
magnetic field H, in the limit where the spin-spin exchange coupling becomes infinitely 
strong, so that the only spin configurations contributing to this partition function are those 
for which spins on adjacent vertices are different There has been continuing interest in 

the Potts model and chromatic and Tutte polynomials for many years; reviews of the Potts 
model include jsjj-jsl and reviews of chromatic and Tutte polynomials include 0]- 15]. 

There are a number of motivations for this study. One is the intrinsic mathematical inter- 
est of these two new kinds of graph coloring problems. A second stems from the equivalence 
to the statistical mechanics of the Potts antiferromagnet in a magnetic field. A third is the 
fact that these weighted graph coloring problems have physical applications. For example, 
the weighted graph coloring problem with < w < 1 (i.e., the DFCP) describes, among 
other things, the assignment of frequencies to commercial radio broadcasting or wireless 
mobile communication transmitters in an area such that (i) adjacent stations must use dif- 
ferent frequencies to avoid interference and (ii) stations prefer to avoid transmitting on one 
particular frequency, e.g., because it is used for data-taking by a nearby radio astronomy 



antenna [16 1. The weighted graph coloring problem with w > 1 (i.e., the FCP) describes this 
frequency assignment process with a preference for one of the q frequencies, e.g., because it 
is most free of interference. 

We remark on some special cases of Ph(G, q, w). The case w — 1 is equivalent to the usual 
(unweighted) chromatic polynomial, P(G,q), counting the number of proper g-colorings of 
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the vertices of G: 



Ph(G,q,l) = P(G,q) . 



(1.1) 



The chromatic number of G, denoted x{G), is the minimal number of colors for which one 
can carry out a proper g-coloring of the vertices of G. For w = 0, one is complete forbidden 
from assigning the disfavored color to any of the vertices, so that the problem reduces to 
that of a proper coloring of the vertices of G with q — 1 colors without any weighting among 
these q — 1 colors, which is thus described by the usual unweighted chromatic polynomial 



Thus, the DFCP, described by Ph(G, q, w) with < w < 1, may be regarded as interpolating 
between P(G, q) and P(G, q — 1). In the FCP, as w increases above 1 to large positive values, 
the favored weighted of one color is increasingly in conflict with the strict constraint that 
no two adjacent vertices have the same color. Hence, the FCP involves frustration in the 
technical sense of statistical mechanics, i.e. mutually conflicting tendencies built into the 
system. 



Consider a graph G = (V, E), defined by its set of vertices V and edges (= bonds) E. A 
spanning subgraph G' C G is defined as the subgraph containing the same set of vertices V 
and a subset of the edges; G' = (V, E') with E' C E. For a graph G we denote the number 
of vertices, edges, and connected components as n(G), e(G), and k(G), respectively. Where 
no confusion can result, we shall often abbreviate n(G) as simply n. We further denote 
the connected subgraphs of a spanning subgraph G' as G[, i = l,..,k(G'). To obtain an 
expression for Ph(G,q,w), we make use of the fact that it is a special case of the partition 
function for the g-state Potts model in an external magnetic field in the limit of infinitely 
strong antiferromagnetic spin-spin coupling. In thermal equilibrium at temperature T, the 
general Potts model partition function is given by 



P(G,q-l): 



Ph(G,q,0) = P(G,q-l) . 



(1.2) 



II. DEFINITIONS AND SOME BASIC PROPERTIES 



A. Relation with Potts Model in an External Magnetic Field 




(2.1) 



with the Hamiltonian 




(2.2) 
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where i, j, £ label vertices (sites) in G, Oi = l,...,q are classical spin variables on these 
vertices, (3 = {k B T)~ l , and (ij) denote pairs of adjacent vertices. Without loss of generality, 
we have taken the magnetic field H to single out the spin value cr» = 1. Let us introduce the 
notation 

K = (3J, h = (3H, y = e K , v = y-l, w = e h . (2.3) 

Thus, the physical ranges of v are v > for the Potts ferromagnet, and — 1 < v < for 
Potts antiferromagnet. The weighted chromatic polynomial is then obtained by choosing the 
antiferromagnetic sign of the spin-spin coupling, J < and taking K — > — oo while keeping 
h = f3H fixed. Since K = j3J, the limit K — > — oo results if one takes J — ► — oo while holding 
T and if fixed and finite. Alternatively, the limit K — > — oo can be obtained by taking the 
zero-temperature limit T — > 0, i.e., /3 — > oo, with J fixed and finite and H — > so as to keep 
h = (3H fixed and finite. The limit .fT — ► — oo guarantees that no two adjacent spins have 
the same value, or, in the coloring context, no two vertices have the same color. One sees 
that in this statistical mechanics context, it is the external magnetic field that produces the 
weighting that favors or disfavors a given value for the spins ov Positive H gives a weighting 
that favors spin configurations in which spins have a particular value, say 1, or equivalently, 
vertex colorings with this value of the color assignment, while negative H disfavors such 
configurations. For positive and negative H, the corresponding ranges of w are w > 1 and 
< w < 1, respectively. 

The partition function Z can be written in a manner that does not make explicit reference 
to the spins <Ji or the summation over spin configurations, but instead as a sum of terms 
arising from spanning subgraphs G' E G. The formula, obtained by F. Y. Wu, is [l^] 

^ fe(G') 

Z(G,q,v,w)= v e{ - G>) |~| [q — 1 + w n ^ . (2.4) 

G'CG i=l 

This can be understood by writing Eqs. (12. ip with (12.21) as 

z =E[n( i+ ^)][n e "M • ( 2 - 5 ) 

{o-n} (ij) I 

If h = 0, then each edge of a particular G' gives a contribution of v and represents a spin 
configuration in which the spins on the ends of this edge have the same value. The spins in 
each component of G' are connected by edges, so they all have the same value, and there 
are q possibilities for this value. In this case, from Eq. (12. 5p one sees that the resultant 
term in summand of (12.41) is simply If h ^ 0, all of the spins in each connected 

component G[ of G' have either the value a,; = 1 or all of these spins have one of the other 
q — 1 values. If they all have the value 1, then each vertex in this G\ gives a contribution of 
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w, so from G\ one gets the contribution w n ( G *\ while if they all have one of the other q — 1 
values, the contribution is 1. In general, therefore, the contribution of the component G\ in 
G' is (q — 1 + w n( - G '^). Taking account of all of the k(G') components in each G' gives the 
factor fJiS^ \ Q ~~ 1 + w n ( G '^j, which is then summed over all G 1 C G. The Wu formula 
( 12. 4p is a generalization of the Fortuin-Kasteleyn formula for the zero-field model 181 ] ■ The 
original definition of the Potts model, (12. ip and (12. 2p . requires q to be in the set of positive 
integers N+. This restriction is removed by Eq. (12 ,4p . Eq. (I2.4p shows that Z is a polynomial 
in the variables q, v, and w, hence our notation Z(G, q, v, w). 



In the special case of zero external magnetic fie 
reduction to the Fortuin-Kasteleyn cluster formula 



d, H = 0, whence w = 1, one has the 



18| 



Z(G,q,v,l) = V<G ' ] 9* (G ° • ( 2 - 6 ) 



G'CG 



This zero-field Potts model partition function is equivalent to the Tutte polynomial 
T(G,x,y), defined by 

T(G, x,y)=J2(*- l) fc(G,) - fc(G) (y ~ 1) C(G,) , (2.7) 

G'CG 

where c(G') is the number of linearly independent cycles in G', satisfying c(G') = e(G') + 
k(G')-n(G'), and 

x=l + ~. (2.8) 

v 

(We remark that k(G') — k(G) and c(G') are the rank and co-rank of G'.) The equivalence 
is given by the relation 

Z(G, q, v) = {x- l) k ^ (y - 1)"( G ) T(G, x, y) . (2.9) 
In Ref. [l| we defined a generalization of the Tutte polynomial, 

U(G,x,y,w) = {x - l)- k ^ G \y - l)- n ^ ]T (y - 
KG') 



X 
G'CG 



X 

8=1 



Y[ (xy - x - y + w n{G '^) . (2.10) 



This function satisfies U(G,x,y,w) — (x — l)~ k( - G \y — l)~ n ( G ^Z(G, q, v , w) and reduces to 
the Tutte polynomial if w = 1: U(G,x,y, 1) = T(G,x,y). 

The .fT — > — oo limit that yields the weighted chromatic polynomial is equivalent to 
v = —1, so 

Ph(G,q,w) = Z(G,q, -l,w) . (2.11) 
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Hence, a constructive formula for Ph{G, q, w) is 

KG') 

Ph{G,q,w) = ^(-l) e(G ' } Yl (q - 1 + w n W) . (2.12) 

G'CG i=l 

For the special case h = 0, i.e., w = 1, one thus has the result of Eq. (11.11) . The limit 
h — > — oo, i.e., w — > 0, effectively removes one of the possible values of the dynamical 
variables c^, or equivalently, in the spanning subgraph formula (12.41) . one of the values of q, 
so 

Z(G,q,v,0) = Z(G,q- l,v,l) . (2.13) 

The special case of this for v — — 1 is Eq. (11.21) . It follows that each of the zeros of 
Ph(G, q, 1) = P(G, q) in the complex q plane shifts to the right by one unit as one replaces 
the value w = 1 by w = 0. In the limit as n — > oo, the accumulation set of the zeros, B q also 
is replaced by its identical image shifted to the right in the q plane as one replaces w = 1 by 
w = 0. 



B. Results for Graphs with Loops, Multiple Edges, and Multiple Components 

If G has any loop, defined as an edge that connects a vertex to itself, then a proper 
g-coloring is impossible. This is because such a g-coloring requires that any two adjacent 
vertices have different colors, but since the vertices connected by an edge are adjacent, the 
presence of a loop in G means that a vertex is adjacent to itself. Thus, 

Ph(G,q,w) = if G contains a loop . (2-14) 

Hence, to avoid having Ph(G, q, v) vanish trivially, we shall restrict our analysis in this paper 
to loopless graphs G. Accordingly, in the text below, where G = (V, E) is characterized as 
having a non-empty edge set E ^ 0, it is understood that E does not contain any loops. 

Another basic property of a chromatic polynomial is that as long as two vertices are joined 
by an edge, adding more edges connecting these vertices does not change the chromatic 
polynomial. This is clear from the fact that the chromatic polynomial counts the number 
of proper g-colorings of the vertices of G, and the relevant condition - that two adjacent 
vertices must have different colors - is the same regardless of whether one or more than one 
edges join these vertices. Let us define an operation of "reduction of multiple edge(s)" in 
G, denoted Re(G), as follows: if two vertices are joined by a multiple edge, then delete all 
but one of these edges, and carry out this reduction on all edges, so that the resultant graph 
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Re{G) has only single edges. Then if G is a graph that contains one or more multiple edges 
joining some set(s) of vertices, 

P(G,q) = P{R E {G),q) . (2.15) 

Since the same proper g-coloring condition holds for the weighted chromatic polynomial, we 
have 

Ph{G, q, w) = Ph(R E (G), q, w) . (2.16) 

Moreover, if G consists of two disjoint parts, G\ and G 2) then Ph(G,q,w) is simply the 
product Ph(G,q,w) = Ph(Gi,q,w)Ph(G 2 ,q,w). Hence, without loss of generality, we will 
generally restrict to connected graphs G. 



C. General Structural Properties of Ph(G,q,w) 

Here we prove some general structural properties of Ph(G, q, w) that hold for an abitrary 
graph G. As discussed above, to avoid having Ph(G,q,w) vanish trivially, we take G to 
be loopless, and without loss of generality, we assume that G is connected. We first apply 
the proper g-coloring condition to analyze properties of Ph(G,q,w) for q e N + . Since this 
proper g-coloring condition cannot be met for integer q — 1, x{G) — 1, it follows that 

x(G)-i 

Ph(G, q, w) contains a factor |~J (q — j) ■ (2-17) 

i=i 

Provided that G = (V, E) has at least one edge, i.e., E ^ 0, the proper g-coloring condition 
cannot be satisfied if q — 1. Hence, a corollary of Eq. (I2.17P is 

If E ^ 0, then Ph(G, q, w) contains a factor (q — 1) . (2-18) 

We can show that if w — 0, then the factor (q — 1) is present even if G does not contain 
any edge. Using our previous result that Z(G,q, v,0) = Z(G,q — l,v, 1) and the fact that 
Z(G, q, v, 1) has a factor of q, we obtain the result that Z(G, q, v, 0) contains the factor q — 1 
and hence Ph(G, q, 0) contains a factor (q— 1). More generally, since Ph(G, q, 0) = P(G, q—1) 
and P{G, q — 1) vanishes for integer q = 1, x(G), it follows that 

X(G) 

Ph(G, q, 0) contains a factor J^J (q — j) . (2.19) 

j'=i 
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Substituting q = in (12.41) and using the factorization 

w n (°i) - 1 = ( w - 1) J2 w e (2-20) 

e=o 

proves that 

Z(G, 0, v, w) contains a factor of (w — 1) . (2-21) 
Setting v = — 1, we thus deduce that [l| 

Ph(G,0,w) contains a factor of (w — 1) . (2.22) 

It is convenient to define the notation 

g = g-l, w=w-l. (2.23) 

From Eq. (12 .4p . it follows that we can write Z(G,q,v,w) in several equivalent ways: 

n(G) e(G) n(G) e(G) 

z(G,q,v,w) = ^2 E a r,s,tq r v s w t = E E Ks,tQ r y s w t 

r,t=0 s=0 r,t=Q s=0 

n(G) e(G) n(G) e(G) 

= E E c r,«.*? W = EE dr^ttfv'w* , (2.24) 

r,t=0 s=0 r,t=0 s=0 

where a nSj i, & r , s ,t, <V S ,*, and d nSjt are integers. Some a r>St t and 6 riSji can be negative, but, 
as we showed in the nonzero c rjSj t and c? r)Sjt are positive. From these equations, one 
infers corresponding ones for Ph(G,q,w) by setting v = —1, i.e., y = 0. Note that in 
the polynomial Z(G,q,v,w) = ^"l=o Ss=o brstQ r y s w\ clearly only the terms with s = 
contribute to Ph(G,q,w). 

From Eq. (12 .4p . it is evident that the term in Ph(G,q,w) of maximal degree in q, or 
equivalently, in q, arises from the contribution of the spanning subgraph G' with no edges, 
for which k[G') = n(G). This term is (with n = n(G)) 

(q + w) n . (2.25) 

It follows that 

a n,0,0 = ^n,0,0 = c n,0,0 = <^n,0,0 = 1 • (2.26) 

There are e(G) spanning subgraphs G' with one edge, since there are e(G) ways of choosing 
this edge. Hence (with our restriction to loopless G), the contribution of these G' in (12.41) is 

e(G)v(q + w 2 ){q + w) n ~ 2 . (2.27) 
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Expanding the terms in Eqs. (12.251) and (12.271) in powers of q and w, we find that the term 
in Z(G, q, v, w) of degree n — 1 in q is 

(e(G)v + to) q n - x . (2.28) 

Similarly, expanding the terms in Eqs. (12. 25ft and (12.271) in powers of q and w, we find that 
the term in Z(G, q, v, w) of degree n — 1 in q is 

(e{G)v + n{w - 1)) q^ 1 . (2.29) 

For our analysis below and for comparisons with chromatic polynomials, it will be con- 
venient to write Ph(G,q,w) as a polynomial in q with w-dependent coefficients, which we 
denote ac^w): 

n 

Ph(G, q,w) = J2 aG,n-jW<l n ~ j • ( 2 -30) 

3=0 

From our discussion above, we have 

a G , n = 1 , (2.31) 

and, using also Eq. (I2.16p . 

a G ,n-i = - (e(R E (G)) + n(l - w)) . (2.32) 

Moreover, from Eq. (12.221) . it follows that the q° term in Ph(G,q,w) contains a factor of 
(w - 1), i.e., 

«g,o contains a factor of (w — 1) . (2.33) 

It is also useful to express Ph(G, q, w) as a polynomial in w with g-dependent coefficients, 
which we denote @G,z{q) (there should not be confusion with /3 = l/(fceT)): 

d w (G) 

Ph(G,q,w)= Yl PgM^ , ( 2 -34) 

3=0 

where d w {G) = deg w (Ph(G, q, w)) is the (maximal) degree of Ph(G, q, w) in w. This degree, 
d w (G) is a G-dependent number less than n. To understand this, we recall that the maximum 
degree of Z{G 1 q ) v ) w) in w is n. This term is y e ^w n = (v + l) e ^w n and corresponds to 
all of the vertices having the same color, 1. However, the possibility that all of the vertices 
have the same color, and, indeed, the possibility that any adjacent vertices have the same 
color, are excluded for Ph(G,q,w), as is evident from the fact that the coefficient of w n 
vanishes for v — — 1. Hence, d w (G) < n. We shall give this degree below for various families 
of graphs. Since all of the nontrivial graphs G = (V, E) that we shall consider have at least 
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one edge, i.e., E ^ 0, Eq. ( 12. 18ft shows that for these, Ph(G, q, w) has the factor (q — 1). In 
analyzing zeros of Ph(G, q, w) it will be convenient to separate this factor out, and we thus 
define, for graphs containing at least one edge, 

PgAq) = (<i-i)PgM > ( 2 - 35 ) 

so that 

d w (G) 

If E ^ 0, then P/i(G, g, w) = (g - 1) ^ PgAq)^ ■ ( 2 - 36 ) 

where f3c,j{q) are polynomials in g. From Eq. (I2.34p and Eq. (11.21) . we obtain the relation 

Ph(G, g, 0) = p G M = P(G, g - 1) . (2.37) 
Combining this with Eq. (12.191) . we have the result that 

X(G) 

Pg,o contains a factor J^J (g — j) . (2.38) 

j=i 

Now, Ph(G,q,l) = T.'jZ^ Poj, but also Ph(G,q,l) = P(G,q), so, using the fact that 
P(G, q) = for integer q = 0, x{G) — 1, we derive the following factorization property for 
the sum of the /3gj coefficients: 

d w (G) x(G)-l 

contains a factor j~J (g — j) . (2.39) 

3=1 3=0 

D. Absence of Deletion-Contraction Relation 

For a graph G, we denote the graph obtained by deleting an edge e £ E as G — e and 
the graph obtained by identifying the two vertices connected by this edge e as G/e. The 
chromatic polynomial satisfies the deletion-contraction relation P(G, q, v) = P(G — e, q) — 
P(G/e,q). In contrast, for w not equal to 1 or 0, the polynomial Ph(G,q,w) does not, in 
general, satisfy this deletion-contraction relation. It is of interest to examine the quantities 
that measure the deviation from such a relation, namely 

APh(G, e, q, w) = Ph(G, q, w) - \ph(G - e, q, w) - Ph(G/e } q, w)] . (2.40) 

We know that APh(G,e,q,w) contains a factor w(w — 1) since for w — 1 and w — 0, 
Ph(G,q,w) is equal, respectively, to P(G,q) and P(G,q — 1), both of which do satisfy the 
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deletion-contraction relation. Furthermore, because of Eq. (12.181) . if G, G — e, and G/e 
contain at least one edge, then APh(G, e, q, w) contains the factor (q — 1). 

As an illustration, using our explicit calculations given below for n-vertex line graphs L n 
and circuit graphs C n , we find the following results. For the first two graphs, L 3 and C3, the 
deletion and contraction on any edge gives the same result, so we need not specify which 
edge is involved. We find, for any edge e, 

APh(L 3 , e, q, w) = APh(C 3 , e, q, w) = -w{w - l)(q - 1) (2.41) 

and 

APh(C 4 , e, q, w) = -w(w -l){q- l)(q - 2) . (2.42) 
For L 4 , denoting e oute r as either of the two outer edges and e m id as the middle edge, we find 

APh(L A , e mid , q, w) = -w{w - l)(g - l) 2 (2.43) 

and 

APh(L 4 ,e outer ,q,w) = -w(w-l)(q-l)(q + w-2) . (2.44) 

It is straightforward to calculate similar differences APh(G,e,q,w) for graphs with more 
vertices and edges, but these are sufficient to illustrate the absence of a usual deletion- 
contraction relation for the weighted chromatic polynomial. 

E. T, P, U, and Ph Equivalence Classes 

An important property of the weighted chromatic polynomial Ph(G, q, w) is the fact that 
it can distinguish between certain graphs that yield the same chromatic polynomial P(G, q). 
More generally, an important property of the partition function of the Potts model in a 
nonzero external magnetic field, Z(G, q, v, w), or equivalently, the function U(G, x, y, w) that 
we defined in Ref. |l|], is that Z(G, g, v , w) and U(G, x, y, w) distinguish between graphs that 
yield the same zero-field Potts model partition function, Z(G,q,v, 1) or equivalently, Tutte 
polynomial T(G, x, y). We begin with some definitions. Two graphs G and H are defined as 
(i) Tutte-equivalent (T-equivalent) if they have the same Tutte polynomial, or equivalently, 
the same zero-field Potts model partition function, Z(G, q, v, 1); (ii) ^/-equivalent if they 
have the same Z(G,q,v,w); (iii) chromatically equivalent (P-equivalent) if they have the 
same chromatic polynomial, P(G, q), and (iv) P/i-equivalent if they have the same weighted 
chromatic polynomial, Ph(G,q,w). 

Let us give some examples. Recall the definition that a tree graph is a connected graph 
that contains no circuits (cycles). The set of tree graphs {G tr ee,n} with a given number, 
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n, of vertices, forms a Tutte equivalence class, with T(G tr ee,n,x,y) = x n_1 , or equivalently, 
Z(G, q, v, 1) = q(q + v) n ~ 1 . However, the Potts partition function in a field, Z(G, q, v, w), or 
equivalently, the function U (G, x, y, w) is able to distinguish between different tree graphs in 
a Tutte-equivalence class. For instance, consider the n = 4 line graph L 4 and star graph £4 
(the graph with one central vertex connected to three outer vertices by corresponding edges). 
These have the same Tutte polynomial T(L 4 , x, y) = T(S4,x,y) = x 3 , or equivalently, the 
same zero-field Potts partition function Z(L 4 , q, v, 1) = Z(Sa, q, v, 1) = q(q + v) 3 , but the full 
Potts partition functions, Z{L^q,v,vo) and Z(S 4 , q, v , w) are different (see Eqs. (I3.6P and 
(I3.16P below). Similarly, L 4 and S4 are chromatically equivalent, with P(L 4 , q) = P(>S 4 , q) = 
q(q — l) 3 as a special case of the result P(G tree>n , q) = q(q — l)™^ 1 for any tree graph with 
n vertices, G tre e,n- However, from our calculations given below in (13.151) and (13.191) . we find 
that P/i(L 4 , q, w) and P/i(S' 4 , q, w) are different. We reach the same conclusion for all of the 
tree graphs that we have studied, i.e., although the set of tree graphs with a given number, 
n, of vertices, forms a chromatic equivalence class, these graphs have different weighted 
chromatic polynomials. We will illustrate this below for n = 5 and n = 6. 

A second set of examples involves graphs with multiple edges. Let us assume that G 
contains one or more multiple edges joining pair(s) of vertices. Such graphs are not Tutte- 
equivalent, but, as noted above, are chromatically equivalent. Because the same proper 
g-coloring condition also holds for weighted chromatic polynomials, these graphs are also in 
the same P/i-equivalence class, as was stated in Eq. (I2.16p . A simple example is provided by 
the line and cirtuit graphs with n = 2 vertices, L2 and C2, the latter of which has a double 
edge connecting the two vertices. One has 

Z(L 2 , q,v,w) = {q-l + wf + v(q - 1 + w 2 ) (2.45) 

and 

Z(C 2 , q,v,w) = (q-l + wf + v(v + 2)(q-l+ w 2 ) , (2.46) 

so that 

Z{C 2 , q, v, w) - Z{L 2l q, v, w) =v(v + l)(q-l + w 2 ) . (2.47) 

The fact that the difference in Eq. (12.471) vanishes for v = —1, i.e., that Ph(L 2 ,q,w) = 
Ph(C 2 , q, w), is a special case of the general result (I2.16p . 

Because of the above-mentioned result that all n-vertex tree graphs are chromatically 
equivalent, in conjunction with the property that Ph(G, q, w) is a chromatic polynomial for 
w = 1 and w = 0, it follows that the difference between Ph(G, q, w) and Ph(H, q, w) between 
two chromatically equivalent graphs G and H must vanish if w = 1 or w = 0. Since these 
are all polynomials, it thus follows that the difference Ph(G, q, w) — Ph(H, q, w) must have 
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w and to - 1 as factors. Furthermore, if q = 1, then 

Z(G,l,v,w) =y e{G) w n(G) 



(2.48) 



If G has at least one edge, then Z(G, 1, v, w) — if y — 0, i.e., v = —1, Now in order to 
be chromatically equivalent, a necessary condition is that two graphs G and H must have 
the same number of vertices, n(G) = n(H), since the degree in q of P(G,q) is n(G). An 
elementary property of the chromatic polynomial P(G, q), proved by iterative application of 
the deletion-contraction theorem, is that the coefficient of the g n ( G ) _1 term is — e(-R^(G)). 
Therefore, another necessary condition that two graphs G and H be chromatically equivalent 
is that 6(Re(G)) = e(RE(H)). Now recall Eq. (12.181) . according to which if G contains at 
least one edge, then Ph(G, 1, w) = 0. Hence, if G and H are chromatically equivalent, then 
either (i) neither contains any edges, in which case Ph(G,q,w) = Ph(H,q,w) = (q — l + w) n , 
where n = n{G) = n(H), or (ii) if G, and hence H, contains at least one edge, Ph(G, 1, w) = 
Ph(H,l,w) = 0. Hence, if G and H are chromatically equivalent and contain at least 
one edge, then the difference Ph(G,q,w) — Ph(H,q,w) contains the factor (q — 1). These 
results on the factors of Ph(G,q,w) — Ph(H,q,w) for chromatically equivalent graphs will 
be evident in our explicit calculations to be presented below. 

F. On the Weighted Face Coloring Problem for Planar Graphs 

Let us consider a planar graph G = (V, E). We recall that the dual of this graph, G*, is the 
graph obtained from G by associating a vertex of G* with each face of G and connecting these 
vertices of G* with edges that cross each edge of G. There is thus a 1-1 isomorphism between 
the vertices, edges, and faces of G and the faces, edges, and vertices of G*, respectively. A 
proper q coloring of the faces of G* is a coloring of these faces with q colors subject to the 
constraint that no two faces that are adjacent across the same edge have the same color. The 
(usual, unweighted) chromatic polynomial P(G, q) satisfies a duality property, namely that 
P(G, q) counts not just the proper q colorings of the vertices of G, but also, and equivalently, 
the proper q colorings of the faces of G*. By the same duality property, for a planar graph 
G, our weighted chromatic polynomial Ph(G, q, w) describes not just the weighted proper q 
colorings of the vertices of G but also, and equivalently, the weighted proper q colorings of 
the faces of G*. 
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III. CALCULATIONS OF Ph(G,q,w) FOR SOME FAMILIES OF GRAPHS 



In this section we give some illustrative explicit calculations of Ph(G, q, w) for various 
families of graphs. Although we generally consider connected graphs, we note that for the 
graph N n consisting of n vertices with no edges, 

Z{N n , q, v, w) = Ph(N n , q 7 w) = {q-l + w) n . (3.1) 

We recall that a tree graph is defined as a connected graph with no circuits. In the following 
text and in Appendix B we present results for the weighted chromatic polynomials of n- vertex 
tree graphs with n up to 6. 



A. Line Graph L n 

The line graph L n is the graph consisting of n vertices with each vertex connected to the 
next one by one edge. One may picture this graph as forming a line, whence the name. For 
n > 2, the chromatic number is x(L n ) = 2. In [2| we gave a general formula for Z(L n , q, v, w), 
and the special case v — —I determines Ph(L n , q,w). Let us define 

, q + v — 1 w 
Tz,i,o =\ | (3.2) 

q — 1 w(v + 1) 



(3.3) 
(3.4) 




and 



Then 



*i = ( J ) • (3-5) 



Z(L n , q, v, w) = u\ H lt0 (T z , lfi ) n ~ l 8l (3.6) 

and Ph(L n , q, w) = Z(L n , q, —1, to). Since e(L n ) — n — 1, we can apply Eq. (12.321) to deduce 
that 

axn,n-i = 1 + n{w - 2) . (3.7) 
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From our general formula (I3.6p . evaluated at v = — 1 to obtain Ph(L n , q,w), we can 
derive some other corollaries concerning coefficients of Ph(L n ,q,w). The maximal degree of 
Ph(L n , q, w) in w is 

deg w (Ph(L n ,q,w))= [—J , (3.8) 

where here [v] denotes the largest integer less than or equal to v e R. This contrasts with 
the fact that the highest power of w in Z(L n , q, v, w) for v ^ — 1 is n. The reason for this is 
that spin configurations that would yield terms of degrees less than or equal to n and greater 
than the maximum in Eq. (13. 8p are forbidden by the proper g-coloring constraint. If n is 
odd, say n = 2m + 1, the coefficient of the term in Ph(L 2m +i, q, w) of highest degree in w, 
namely the coefficient of the term w^ n+1 ^ 2 = w m+1 , is 



0l 



L2m+l,m+l 



(3.9) 



If n is even, say n = 2m, the coefficient of the term in P/z(L 2m , q, w) of highest degree in w, 

(3.10) 



namely the coefficient of the term w n l 2 = w m , is 

f3 L2m , m = (q-l) m - 1 ((m + l)q- 
The coefficient of the w° term in Ph(L n , q, w) is 

P Lnfi = (q-l)(q-2r~ 1 



2m) . 



(3.11) 



We proceed to give some explicit results for Ph(L n ,q,w) for various values of n. The 
case Li = Ni is already covered by Eq. (I3.ip . For the next few cases we list the explicit 
polynomials below, both in factored form and in the form of Eq. (12.301) : 

Ph(L u q,w) =q-l + w (3.12) 



Ph(L 2 ,q,w) = (q-l) q + 2(w-l) 



q 2 - (3 - 2w)q + 2(1 - w) (3.13) 



Ph(L 3 ,q,w) = (q-l) q 2 + (3w-4)q+(w-l)(w-4) 



q 3 - (5 - 3w)q 2 + (w 2 - 8w + 8)q - (w - l)(w - 4) (3.14) 



and 



Ph(L 4 ,q,w) = (q-l)(q + w- 2) q 2 + (3tw - A)q - 4(w - 1) 
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= g 4 - (7 - 4w)q 3 + 3(w 2 - 6w + 6)q 2 - (7w 2 - 26w + 20)q 



+ 4(u;-l)(tw-2) . (3.15) 
Results for tree graphs with n = 5 and n = 6 vertices are given in Appendix B. 

B. Star Graphs S n 

A star graph S n consists of one central vertex with degree n — 1 connected by edges 
with n — 1 outer vertices, each of which has degree 1. For n > 2, the chromatic number is 
x(S„) = 2. We have derived the following general formula for Z(S n ,q,v,w): 

n-l / _ i \ 

Z(S n ,q,v,w) = J2[ n ■ )^{q + w^ 1 ){q + w) n - 1 -\ (3.16) 

where q was given in Eq. (I2.23p . Evaluating this for v — — 1 yields Ph(S n , q, w). The term 
in Ph(S n ,q,w) of maximal degree in w corresponds to a configuration in which all of the 
outer vertices are assigned the color 1 and the central vertex of the star graph is assigned 
any of the other q — 1 colors. For n > 3 where the star graphs are nondegenerate, this term 
is thus (q — l)u> n_1 , so that, in particular, 

deg w (Ph(S n ,q,w)) = n-l . (3.17) 

(The graph S2 is degenerate in the sense that it has no central vertex but instead coincides 
with L 2 .) The graph S3 is nondegenerate, and coincides with L 3 . For n = 2, the term 
in Ph(S2,q,w) of maximal degree in w, namely the coefficient of the term w, is 2(q — 1). 
For n > 3, the coefficient of the term in Ph(S n ,q,w) of maximal degree in w, namely the 
coefficient of the term w n ~ l , is (q — 1). This is easily understood since it corresponds to the 
assignment of the color 1 to each of the n — 1 outer vertices of the star graph S n , which 
allows any of the remaining (q — 1) colors to be assigned to the central vertex of this graph. 
Because the number of edges of the star graph is e(S n ) — n — 1, it follows that 

= 1 + n(w - 2) . (3.18) 

This coefficient is equal to aL„, n -i- 

As an explicit example, for the graph S4 we calculate 

Ph(S 4 , q, w) = (q- 1) \q 3 + 2(2w - 3)q 2 + (3w 2 - Uw + 12) q + {w - l)(w 2 - 5w + 8) 
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= g 4 - (7 - 4w)q 3 + 3(w 2 - 6w + 6)q 2 - (-w 3 + 9w 2 - 27w + 20)q 

+ (1 -w)(w 2 -5w + 8) . (3.19) 
Results for S n with n = 5 and n = 6 are given in Appendix B. 



C. Distinguishing Between Some Chromatically Equivalent Graphs 



Using the results given in the text and Appendix B for tree graphs with up to six vertices, 
we now analyze the differences between the weighted chromatic polynomials for tree graphs 
that are chromatically equivalent. There are two tree graphs with n = 4 vertices, namely, 
L 4 and S4. In chemical nomenclature, L 4 and S4 correspond to the carbon atom backbones 
of the alkanes n-butane and isobutane (i.e., 2-methylpropane) . From Eqs. (13.151) and (13.191) 
we find 

Ph(Si, q, w) - Ph{L 4 , q, w) = (q- l)w(w - l) 2 . (3.20) 



L 5 



S 5 



FIG. 1: Tree graphs with n = 5. 



There are three tree graphs with n = 5 vertices, as shown in Fig. (TJ namely (i) the line 
graph L 5 , (ii) a graph that we denote Y 5 , which is obtained by starting with the star £4 
and elongating one of the edges by the addition of another vertex and edge, and (hi) the 
star graph S$. We order this list in terms of graphs of increasing maximal vertex degree 
A; one has A = 2, 3, 4 for L 5 , Y 5 , and S5, respectively. In chemical terminology, these 
correspond to the carbon atom backbones of the alkanes (i) n-pentane, (ii) isopentane (i.e., 
2-methylbutane), and (hi) neopentane (i.e., 2,2-dimethylpropane), respectively. From Eqs. 
( IBB . ([BSD , and (EU), we calculate 

Ph{S 5 , q, w) - Ph(L 5 , q, w) = (q- l)w(w - l) 2 (3g + w-S) (3.21) 

Ph(S 5 , q, w) - Ph(Y 5} q, w) = (q- l)w(w - l) 2 (2q + w - 3) (3.22) 
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and thus 

Ph{Y 5 , q, w) - Ph(L 5 , q, w) = {q- l)(g - 2)w(w - l) 2 . (3.23) 

For these graphs, one observes that the differences between chromatically equivalent graphs 
have a double zero at w = 1. We find that this is also true of the differences between weighted 
chromatic polynomials of tree graphs with n = 6 vertices, as discussed in Appendix B. 



D. Complete Graphs K n 

The complete graph K n is the graph with n vertices such that each vertex is connected 
to every other vertex by one edge. The chromatic number is thus x{K n ) — n - One has 
e(K n ) = (™). The simplest two cases coincide with previously discussed graphs, namely the 
single vertex, K\ = L 1( for which we gave Ph(Li, q,w) in Eq. (13.121) . and the n = 2 case, 
for which K 2 = L 2 and Ph(L 2 ,q,w) was given in Eq. (13.131) . For general n > 2 we obtain 
the following theorem: 

n-l 

Ph(K n , q,w)=[Y[{q- j)] (q + n(w - 1)) . (3.24) 

i=i 

Proof: To prove this, we begin by observing that because of the proper g-coloring condition, 
Ph(K n ,q,w) vanishes for all of the integer values q = 1, ...,n — 1 and hence must contain 
the factor YllZi {q ~ j)- The proper g-coloring condition also means that only one vertex at 
most can be assigned the color 1; hence the term in Ph(K n ,q,w) of highest degree in the 
variable w has degree 1. Since the maximal degree of Ph(G,q,w) in the variable q is n(G), 
it must be of the form 

n-l 

[i[(q-3)](aq + bw + c) . (3.25) 

3=1 

From Eq. (12.261) . it follows that a = 1. From Eq. (II. ip we have 

n— 1 

Ph(K n , q, 1) = P(K n , q) = H(q- j) , (3.26) 

3=0 

which implies that b = — c, so the last factor in (13.251) is (q + b{w — 1)). From Eq. (11.21) we 
have 

n 

Ph(K n , q, 0) = P(K n , q -l) = Y[( q - j) , (3.27) 

3=1 

which implies that b = n, so that the additional factor is (q + n{w — 1)). This proves the 
result (T3T24j) . □ 
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A corollary of the theorem of Eq. ( 13.241) is that 

deg w (Ph(K n ,q,w)) = l (3.28) 
and, further, for n > 2, the term in Ph(K n , q, w) of maximal degree in w has coefficient 



ra-1 



(3.29) 



E. Wheel Graphs Wh n 

The wheel graph VF/i n is the graph obtained by joining one central vertex to the n — 1 
vertices of the circuit graph C n . (This is the "join" of K\ with C n -i.) The central vertex 
can be regarded as forming the axle of the wheel, while the n — 1 vertices of the C n _i and 
their edges form the outer rim of the wheel. This is well-defined for n > 3, and in this 
range the chromatic number is x{Wh n ) = 3 if n is odd and x{Wh n ) = 4 if n is even. The 
graph WI13 involves one double edge, while the Wh n graphs with n > 4 have only single 
edges. The first nondegenerate case is Wh 4 , which is the same graph as K 4 . We have given 
the general structure of Z(Wh n+ i,q,v,w) in Ref. [2j, and this determines the structure of 
Ph(Wh n+ x, q, w). Reductions for w — 1 and w = are given in Refs. [3, Eo|. For the 
nondegenerate cases n > 3, the number of edges is e(Wh n ) = 2(n — 1). We calculate 

Ph(Wh n+1 ,q,w) = (q-l)\(\ Wh>+ ) n +(\ Wh> „) n ] +(g-l)(g-3)(-ir 



+ w 



where 



with 



A 



( g _2)»+(g-2)(-l) B 

1 



Wh,± 



-3± v 7 ! 



Wh 



(3.30) 
(3.31) 



A Wh = (g-3) 2 + 4«;(g-2) . (3.32) 

We note that Awh is equal to A\ (given in Eq. (14.231) ) with q replaced by q — 1, so that the 
eigenvalues \wh,± are the same as the eigenvalues A^oj, j = 1,2 (given in Eq. (14.22)) ) with 
q replaced by q — 1: 

^Wh,± — (Al,0j)g^ 9 -1 > (3.33) 

where ± corresponds to j = 1,2, respectively. The A^oj, j = 1,2, enter in Ph(L n ,q,w), 
given above, and Ph(C n ,q,w), given in Eq. (14.211) . From these observations, it follows that 



Ph{Wh n+1 , q, w) = {q- l)Ph(C n , q-l,w) + wP(C n , q - 1) . 



(3.34) 
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This relation makes the reductions of Ph(Wh n+ i,q,w) for w = 1 and w = obvious; using 
Ph(G, q, 1) = P(G, q), one has 

Ph(Wh n+h q,l) = (q-l)P(C n ,q-l)+P(C ni q-l) 

= P(Wh n+1 ,q) = q[(q-2) n + (q-2)(-l) n ] (3.35) 

and 

Ph(Wh n+h q,0) = (q-l)P(C n ,q-2) 

= P(Wh n+1 , q - 1) = (q - 1) [(g - 3)™ + (g - 3)(-l)»] . (3.36) 
Further, from the values of x(Wh n ) for odd and even n, it follows that 

If n is odd Ph(Wh n ,q,w) contains a factor (q — l)(g — 2) . (3.37) 

and 

If n is even Ph(Wh n ,q,w) contains a factor (q — l)(q — 2)(q — 3) . (3.38) 

Although Wh 3 differs from C3 = K 3 in having one double edge, Eq. (I2.16P shows that 
Ph(Wh 3 ,q,w) = Ph(C 3 ,q,w), where Ph(C 3 ,q,w) was given above in Eq. (14.271) . Fur- 
thermore, the graph Wh^ is the same as K4, so Ph(Whi,q,w) = Ph(K 4 ,q,w), where 
Ph(K n ,q,w) was given above in Eq. (13.241) . 

Since the number of edges in the wheel graph e(Wh n+ i) = 2n, we can apply Eq. (12.321) 
to deduce that 

o£\vh n+1 ,n = - (3n + 1 - (n + l)w) . (3.39) 

For the following we again assume that n > 3 so that the Wh n graph is well-defined. The 
highest power of w in Ph(Wh n , q, w) is 

"77 — 1 ~ 

deg w (Ph(Wh n ,q,w)) = [—J . (3.40) 

If n is even, say n = 2m, then the coefficient of the term in Ph(Wh n ,q,w) of maximal 
degree, namely the coefficient of the term w m ~ l , is 

= (2m - l)(g - l)(g - 2) m -\q - 3) . (3.41) 

If n is odd, say n = 2m + 1, then the coefficient of the term in Ph(Wh n ,q,w) of maximal 
degree, namely the coefficient of the term w m , is 

Pwk 2m+1 ,m = 2(q-l)(q-2r . (3.42) 
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As an illustration, we display Ph(WIi5, q, w) below: 
Ph(Wh 5 , q, w) = (q- l)(q - 2) \q 3 - 5(2 - w)q 2 + (2w 2 - 29w + 34)g - (to - l)(4w - 39) 

= g 5 - (13 - 5w)g 4 + 2{w 2 - 22w + 33)g 3 - (I0w 2 - UOw + 161)g 2 

+ (16w 2 - 187w + 185)? - 2(w - l)(4w - 39) . (3.43) 

IV. P/t(G, g, «;) FOR LATTICE STRIP GRAPHS WITH PERIODIC 
LONGITUDINAL BOUNDARY CONDITIONS 

A. General Structure 

In [l|, 0] we have given a general structural formula for Z(G S , L y x m, BC, q, v, w) on strip 
graphs G s of width L y vertices and length L x , with cyclic (eye.) or Mobius (Mb) boundary 
conditions (BC's). For cyclic strips the special case of this structural formula is 

L y n Ph (L y ,d) 

Ph(G s ,L y xm,cyc. } q } w) = ^ id) E l X G s ,L v ,dj^ W T > ( 4 - X ) 

d=0 3=1 

where m = L x for strips of the square and triangular lattices and m = L x /2 for strips of the 
honeycomb lattice. The coefficients are given by 



cW = Y,(-iy( 2d . j ){q-l) d -i. (4.2) 



j=0 



The first few of these coefficients are = 1, = q — 2, c^ 2 -* = q 2 — 5q + 5, etc. For 
Mobius strips, there is a switching of certain c^'s as specified in general in j^j], using the 



same methods that we employed in [2l|, as specified by Eqs. (2.30)-(2.32) and the v = —1 
special case of Eq. (2.33)) of [2[]). 

The numbers nzh(L y ,d) of A's corresponding to each in the general Potts model 
partition function are reduced for the special case v = —1 of interest here. By coloring 



combinatoric arguments similar to those used in 21] and [2j we determine the nph{L y , d) as 



follows. The numbers riph(L y ,d) are identically zero for d > L y , and 

n Ph (L y ,L y ) = l (4.3) 
n Ph {L y ,L y -l) = 2L y (4.4) 
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n Ph (L y + 1,0)= n Ph (L y , 0) + n Ph (L y , 1) (4.5) 

and, for 1 < d < L y + 1, 

n Ph (L y + 1, d) = n Ph (L y + 1, d - 1) + 2n Ph (L y , d) + n Ph (L y , d + 1) . (4.6) 
The n Ph (L y ,d) satisfy the identity 

^ WnpniLy, d) = P(T Ly , q) = q(q - l)^' 1 . (4.7) 

Indeed, one method of calculating the n P h(L y ,d) is to differentiate this equation L y times. 
One thereby obtain L y +1 linear equations in the L y + 1 unknowns n P h(L y , d), d = 0,1, L y ; 
solving these equations yields the results given above. We note that 

n Ph (L y ,0) = C Ly . 1 + C Ly , (4.8) 

where C n is the Catalan number, 

Cn = — — _ ■ (4.9) 



n + 1 \ n 

(No confusion should result with the use of C n to mean both Catalan number and the cyclic 
graph with n vertices, since the context makes clear which is meant.) We recall that the 
partition function of the zero-field Potts model on cyclic strips of the square lattice (as well 
as other lattices) has the structure 



21 



22] 



L y n z (L y ,d) 

Z(G s ,L y xm,cyc.,q,v) = J2c {d) E t^W^F , (4-10) 

d=0 3=1 

where the coefficients are given by Eq. (14. 2 p with q — > q + 1 and 

for < d < L y . For this h = case, the total number of distinct eigenvalues that enter in 
Eq. (14.101) for a lattice strip, i.e., 

Ly 

N ZtLy = J2^z(L y ,d) , (4.12) 

is 

Nz, Ly = for h = 0. (4.13) 
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TABLE I: Table of numbers nph(L y , d) and their sums, Nph,L y f° r strips of the lattice A (square, triangular, 
or honeycomb). Blank entries are zero. See text for further discussion. 



L„ \ d 
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7 


8 


Nph T 


1 


9 
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8 
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7 
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26 


4 


19 


37 


25 


8 


1 










90 


5 


56 


118 


95 


42 


10 


1 








322 


6 


174 


387 


350 


189 


63 


12 


1 






1176 


7 


561 


1298 


1276 


791 


327 


88 


14 


1 




4356 


8 


1859 


4433 


4641 


3185 


1533 


517 


117 


16 


1 


16302 



We find an interesting relation connecting the numbers nph(Ly, d) with the corresponding 
numbers riz{L y ,d) for the zero-field Potts model, namely, for L y > 2, 

n Ph (L y , d) = n z {L y , d) + n z {L y - 1, d) . (4-14) 

From our determination of the n Ph (L y , d), we next calculate the total number 

Ly 

N PhtLy = Yj n Ph{L y ,d) . (4.15) 

d=0 

We find 

2L y \ (2(L y — 1 



^UJ + U-iJ- (416) 

From the relation (14.141) it follows that the total number Nph,L y satisfies the relation, for 
Ly > 2, 

N Ph , Ly = N ZiLy + N z>Ly ^ , (4.17) 

as is evident in Eq. (14 .161) . We list the n Ph (L y , d) and N Ph>Ly for 1 < L y < 8 in Table fl] For 
purposes of comparison, we include tables of n P (L y , d), N P)Ly1 and n z (L y , d) for both h = 
and h (from |2j) in Appendix A. 
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Let us denote Nz,l v for /i 7^ as N Z h,L y for notational clarity, to distinguish this from 
Nz(L y , d) for h = 0. For the Potts model in a nonzero field, we have found p], Q 

^ = E(^)(j-) for MO. (4.18) 

Concerning the relative sizes of Np t L y , Nph t L y , ^z,L y and Nzh,L y , we have, for = 1, 
^Vp,i = Nz,i = 2 < iVpft,! = iV^i = 3 and the inequality 

N P , Ly < N z , Ly < N PhtLy < N zh , Ly for L y > 2 . (4.19) 

For example, this set of four numbers is (4,6,8,11) and (10,20,26,45) for L y = 2 and L y = 3, 
respectively. 

For large strip width L y , Np^L has the same general asymptotic behavior as Nz t L y '- 

N Ph ,L y ~ const, x L y ~ 1/2 4 Ly as L y -> oo . (4.20) 



B. Circuit Graphs C n 

The circuit graph C n , or equivalently, the ID lattice with periodic boundary conditions, 
has chromatic number xiPn) = 2 if n is even and x(C n ) = 3 if n > 3 is odd. (The case 
n = 1 is a single vertex with a loop, for which there is no proper g-coloring, so Ph(C\, q, w) 
vanishes identically.) The polynomial Ph(CL, q, w) for the circuit graph C n , can be obtained 
from the calculations of Z{C n ,q,v,w) 0, Q by setting v = —1. Expressed in our present 
notation, it is 

Ph(C n ,q,w) = (A 1A i)" + (A 1A2 ) n + (q~ 2)(A lti r , (4.21) 

where 

A w = ^[g-2±V^] (4-22) 
where the ± sign corresponds to j = 1,2, 

A 1 = (g-2) 2 + 4(g- l)w (4.23) 

and 

Ai,i = -1 . (4.24) 

From the values of x{C n ) given above, it follows that, in addition to the general factor of 
(q— 1) present for any n, if n is odd, Ph(C n , q, w) contains a factor of (q— 2). Since e(C n ) = n, 
it follows that for n > 3 

«c„,n-i = n(w - 2) . (4.25) 



24 



(For n = 2, C2 has a double edge, so one uses Eq. ( 12. 32ft to obtain ac 2 ,i 
exhibit Ph(C n , q, w) for 2 < n < 5 below: 



-(3-2io).) We 



P/i(C 2 , g, 10) = Ph(L 2 , g, w) = (q- l)[q + 2(w - 1)] = q 2 - (3 - 2w)g + 2(1 - w) (4.26) 



Ph(C 3 ,q,w) = (g-l)(g-2) g + 3(^-1) 



g 3 - 3(2 - w)g 2 + (11 - 9w)q - 6(1 - w) 



(4.27) 



P/i(C 4 , q, w) = (g - 1) q 3 + (Aw - 7)q 2 + (2w 2 - 16w + 17) q - 2(w - l)(w - 7) 



q 4 - 4(2 - w)q 3 + 2(w 2 - lOw + I2)q 2 - (Aw 2 - 32w + 31)g 



+ 2(w -\)(w -1) 



(4.28) 



and 



Ph(C 5 ,q,w) = (q-l)(q-2) q 3 + (5w - 7)q 2 + (5w 2 - 20w + 17) q - 5(w - l)(w - 3) 



g 5 - 5(2 - w)q 4 + 5(w 2 - 7w + 8)q 3 - I0(2w 2 - 9w + 8)q 2 



+ (25w 2 - lOOw + 79)g - 10(w - l)(w - 3) . 



(4.29) 



In the special case of zero-field, h = 0, i.e., w — 1, Ai i0 ,i — q — 1 while Ai ,2 becomes 
equal to A14. Thus, a "transmigration" process occurs in which one of the A's associated 
with the coefficient of degree d = becomes equal to, and hence can be grouped with, a 
A associated with a coefficient with a different degree <i, here d = 1. Hence, one has the 
reduction 

Ph(C n , q, 1) = P(C n , q) = (q~ l) n + (q - l)(-l) n . (4.30) 

For w = 0, Ai >0 ,i = g — 2 while Ai ,2 = 0, so that Ph(C n , q, 0) = Ph(C n , q — 1, 1), in agreement 
with the general relation (11.21) . 

As an application of our result (I2.18P above, it follows that Ph(C n , q, w) contains the factor 
(g — 1). We note some additional factorization properties and special values of Ph(C n , g, w): 



If n is odd, then Ph(C n , g, w) contains the factor (g — 2) 
For the g = 2 case, 



Ph(C n ,2,w) 



i + (-i: 



w 



n/2 



(4.31) 
(4.32) 
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The highest power of w in Ph(C n , q, w) is 



deg w (Ph(C n ,q,w)) 



(4.33) 



where here [v\ denotes the integral part of v. This contrasts with the fact that the highest 
power of w in Z(C n , q, v, w) for v ^ — 1 is n. The reason for this is that spin configurations 
that would yield terms proportional to w p with n/2 < p < n for even n and (n— \) /2 < p < n 
for odd n are forbidden by the proper g-coloring constraint. For n even, say n = 2m, the 
term in Ph(C2m,q,w) of maximal degree in w, namely w m , has coefficient 2(q — l) m . For 
n odd, say n = 2m + 1 with m > 1, the term in Ph(C2 m +i,q,w) of maximal degree in w, 
namely w m , has coefficient (2m + l)(q — l) m (q — 2). 

If and only if w — 1, then Ph(C n , q, 1) = P(C n , q) contains g as a factor. For this zero- 
field case w = 1 we also recall that Ph(C n , q, 1) = P(C n , q) contains q(q — 1) as a factor and, 
furthermore, if n > 3 is odd, then Ph(C n , q) also contains (q — 2) as a factor. 



We denote the cyclic and Mobius strips of the square lattice of width L y = 2 and length 
L x = m as the ladder graph L m and the Mobius ladder graph ML m . For both of these our 
general structure determination above gives nph(2,0) = 3, nph{2, 1) = 4, and nph(2,2) = 1, 
for a total of = 8 terms. The weighted coloring polynomial for L m is 



C. 



L y = 2 Cyclic Strip 



2 



np h (2,d) 




rn 



(4.34) 



where 



A 2 ,o,i = w(l - q) 



(4.35) 




(4.36) 



A 2 = q 4 + 6q 3 w + q 2 w 2 - I0q 3 - 36q 2 w - 2qw 2 + 39q 2 



+ 72qw + w 2 - 70q - 50w + 49 



(4.37) 




(4.38) 



A 3 = q 2 + A(q-l)(w-l) 



(4.39) 
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^2,1, j = 

A, 



and 



q - 4 ± y/A 4 
- Aq(w - 2) + 

1 . 



, J = 3,4 
4(4 - 3w) 



(4.40) 
(4.41) 



A 2 , 2 = 1 • (4.42) 

In Eq. (14 .40 p . j — 3 and j = 4 apply for the + and — sign choices, respectively, and 
similarly for the other equations. Results for Z(L mj q,v,w) are given in Refs. [^J and [^j. 
The weighted chromatic polynomial for the L y = 2 Mobius strip of the square lattice is 
obtained by applying the results of Ref. 
For w — 1, these A's reduce as follows: 



q 



^2,0,1 

^2,0,2 - > 
-^2,0,3 _ 
^2,1,1 _ 
^2,1,2 
^2,1,3 _ 
^2,1,4 



l-q 

-3g + 3 
3-q 
l-q 
1 

3-q 
-> 1 . 



(4.43) 

(4.44) 
(4.45) 
(4.46) 
(4.47) 
(4.48) 
(4.49) 



Thus, a transmigration process of A's occurs here, just as it did for the L y = 1 case; (i) one 
of the three A's in the d = subspace reduces to the single A, q 2 — 3q + 3, in the d = 
subspace for the chromatic polynomial P(L m ,q), while the other two become equal to the 
two A's in the d = 1 subspace of P(L m ,q); (ii) two of the four A's in the d = 1 subspace 
reduce to the two A's, 3 — q and 1 — q, in this subspace for P(L m ,q), while the other two 
become equal to the single A = 1, in the d = 2 subspace for P(L m ,q). Hence, we have the 
reduction 



Ph{L m , q, !) = (!- q) rn + {q 2 - 3q + 3) m + (3 - q) m + c« (1 - q) m + 2 + (3 - q) 



(q 2 -3q + 3) m + c (1) (3 - q) m + (!-?) 



+ c 



(2) 



(4.50) 
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V. SOME PROPERTIES OF THE ZEROS OF Ph(G,q,w) 



A. Zeros of Ph(G, q, w) in q as Functions of w 



Here we discuss the zeros of Ph(G, q, w) in q as a function of w for some illustrative graphs 
G. Since the maximal degree of Ph(G, q, w) in the variable q is n(G), it has this number of 
zeros in the variable q. In contrast, as is evident from our explicit calculations above, the 
maximal degree of Ph(G, q, w) in the variable w depends on details of G. As is true for any 
polynomial, the positions of the zeros of Ph(G, q, w) are continuous functions of q for fixed 
w and continuous functions of w for fixed q. As noted above, for any graph G with at least 
one edge, Ph(G, q, w) contains the factor (q — 1), so it has a fixed zero at q = 1. A general 
statement is that since Ph(G, q, 1) = P(G, q) and Ph(G, q, 0) = P(G, q — 1), it follows that 
each zero of Ph(G, q, w) shifts horizontally to the right by one unit in the complex q plane 
if one replaces w = 1 by w = 0. 

One relevant quantity of interest is the maximal real zero of Ph(G, q, w), which we denote 
q-mrziG). This is related to chromatic number of the graph G, x(G), because, by the definition 
of q m rz{G), Ph(G,q,w) is nonzero for real q > q mrz {G), and it must be positive since for a 
given G and w, if q is sufficiently large, then Ph(G, q, w) is positive. Hence, in addition to 
x{G), which is fixed for a given G, the quantity q mTZ {G) serves as a w-dependent measure of 
the ability to perform a proper vertex coloring of this graph. As a corollary of the discussion 
above, a general result is that q mr z{G) shifts one unit to the right as w decreases from 1 to 
0, as a consequence of (11.11) and (11.21) . 

Let us consider some simple examples. We have 



This increases from to 1 as to decreases from 1 to in the DFCP interval and decreases 
from through negative values as w increases above 1 in the FCP interval. For L 2 , 



This increases from to 2 as w decreases from 1 to in the DFCP interval and decreases 
from through negative values as w increases from 1 in the FCP interval. This example 
also illustrates how the multiplicity of zeros can change as a function of w; for w = 1/2, 
Ph(L 2 , q,l/2) has two coincident zeros at q = 1. 

For L 3 , the situation is more complicated. The expression for Ph(L 3 , q, w) is given in Eq. 
(I3.14p . In addition to the fixed zero at q = 1, Ph(L 3 , q, w) has two other zeros, which occur 



q m rz{.Li) = 1 - w . 



(5.1) 



qmrz(L 2 ) = 2(1 - w) • 



(5.2) 
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at the values 



QL3z,± — 7) 



4 - 3w ± y/w(5w - 4) 



(5.3) 



For w — 1, these reduce to g = 1 and g = for the ± signs. As w decreases from 1, the zero at 
q = increases while the zero at q — 1 decreases. As w decreases through the value w — 4/5, 
these two zeros meet at g = 4/5, and then move off the real axis as a complex-conjugate 
pair as w decreases further in the interval < w < 4/5. The magnitudes of the imaginary 
parts of these complex zeros increase to maximal values as w decreases through the value 
w = 2/5 and then decrease toward zero. As w decreases through the value w — 0, these 
zeros return to the real axis, becoming a double zero at q — 2. Thus, for w — 1, qi3z,+ = 1, 
while for w = 0, g_L3 2j + = 2. However, qi3z,+ does not increase monotonically from 1 to 2 as 
w decreases from 1 to 0; instead, it actually decreases from 1 to 4/5 as w decreases from 1 
to 4/5, while qi3z- increases from to 4/5. As w decreases below 4/5, qL3z,± form a pair 
of complex-conjugate roots, as noted. Hence, for w in the DFCP interval, the fixed zero at 
q — 1 is q mr z(Ls). As w passes through the value w — 0, q rmz (L 3 ) jumps discontinuously 
from the fixed zero at q — 1 to qL3z,+ = Ql3z- = 2. 

As w increases above 1 in the FCP interval, qL3z,+ decreases monotonically from 1, so that 
q rmz (Li) remains the fixed zero at q — 1. This example shows that although individual zeros 
of Ph(G,q,w) in the q plane are continuous functions of w, the maximal real zero q rm z{G) 
of Ph(G, q, w) is a discontinuous function of w. The reason for the discontinuity in q rmz (G) 
is the confluence of two complex-conjugate roots that come together and pinch the real axis 
(at q = 2) as w decreases through w — 0, abruptly producing a new maximal real root (of 
multiplicity 2). Thus, 



but q rm z(L 3 ) = 2 for w = 0. 

In the FCP interval w > 1, the fixed zero at q — 1 remains as q mrz (L 3 ), since 
Ql3z,+ decreases below 1, while qL3z- decreases below 0. Indeed, the zero qL3z,+ has a 
local maximum at w — 1 and decreases monotonically as w increases above 1; qL3z,+ 
passes through as w increases through the value w = 4 and behaves asymptotically like 
QL3z,+ ~ —(1/2) (3 — \fb jtoasiin oo. The zero qi3z,- also decreases monotonically from its 
value of at w — 1 through negative values as w increases above 1, and has the asymptotic 
behavior qi3z- ~ _ (l/2)(3 + V§ )w as w — > oo. 

As u> decreases through negative values, the double zero qL3z,± at g = 2 splits apart again. 
The zero qL3z,+ increases monotonically as w decreases through negative values, and grows 
asymptotically as qL3z,+ ~ (l/2)(3 — V5 )\w\ as w — > — oo. The other zero, gL3 2) - is a 
non-monotonic function of it first decreases below 2, reaching a minimum of qLz3- = 9/5 



lim q rmz {L 3 ) = 1 , 



(5.4) 
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for w = — 1/5 and then increases, passing through the value qhzz- = 2 again as w decreases 
through the value w = — 1, and increasing asymptotically as qizz,- ~ (l/2)(3 + V§ )\w\ 
as w — > — oo. As these examples show, there is somewhat complicated behavior of the 
individual zeros as a function of w for even a very simple graph such as L 3 , and this behavior 
is, understandably, more complicated for larger graphs. 

In this case and others one can avoid the discontinuous behavior of q rm z{G) by restricting 
w to the range w > 0. Doing this, we have examples from these simple graphs that exhibit 
a continuous increase of q rm z{G) as w decreases from 1 to and an example in which q rmz is 
fixed, independent of w in this DFCP interval < w < 1. For w in the FCP interval w > 1 
we find cases where q rm z{G) decreases monotonically as w increases and also a case where 
qrmz(G) is fixed at 1. 

Among the Ph(C n , q, w) polynomials, the case n = 2 is the same as Ph(L,2, q, w) and for 
n = 3, there are two w-independent zeros, at q = 1 and q = 2, while the third occurs at 
q = 3(1 — w) . As w decreases from 1 to 0, this third zero increases from to 3. In the FCP 
range w > 1, this zero decreases from at w = 1 to — oo as iu ^ oo. 

A particularly simple case to discuss is that of complete graphs K n . For these, as is 
evident from Eq. (13.241) . Ph(K n , q, w) has zeros in q at q — 1, 2, .., n — 1, and g = n(l — u>). 

The zeros of Ph(G, q, w) in g as a function of w for fixed w satisfy certain boundedness 
properties [25(. For w — 1 and u> = 0, these specialize to the bound for a (usual, unweighted) 
chromatic polynomial, namely that if q is a zero of P(G,q), then |g| < aA max (G), where 
A max (G) denotes the maximal degree of the vertices in G and a ~ 7.964 from Ref. 25] 
(improved in Ref. |26j). However, these zeros of Ph(G, q, w) in q are unbounded as \w\ — > oo. 



This is already evident in the simplest case of a single vertex, for which Ph(Li,q,w) = 
q — 1 + w, with a zero at q = 1 — w with a magnitude that goes to infinity as \w\ — ► oo. 
Some insight into the lack of boundedness of the zeros of Ph(G, q, w) for arbitrary w can 



be gained by examining the behavior of the factor Hi=i i (Q ~~ 1 + w n ' G ''') in Ph(G,q 
As if — > oo, each of these factors goes to infinity also unless q behaves like 1 — w n ( G '^, 
going to — oo. In principle, one might imagine a cancellation occuring between different 
rii=i*(9 — 1 + w n( - G »- ) )factors for different spanning subgraphs G' C G, via different signs 
of the (— l) e ( G ') factor in Ph(G,q,w), nevertheless, this makes is understandable why, in 
the absence of such cancellation, some zero(s) of Ph(G,q,w) have magnitudes \q\ —>■ oo as 
\w\ — > oo. Note that this behavior cannot simply be attributed to the frustration that occurs 
when w gets large and positive, because it is also true in the unphysical region for w negative. 

Although the positions of the zeros of Ph(G, q, w) in q are continuous functions of w and 
vice versa, this is not true of the asymptotic locus B. Indeed, we shall show below that 
for the n — > oo limit of the circuit graph C n , as w decreases below one, regardless of how 
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small the magnitude of 1 — w is, the part of B q that crosses the real q axis on the left jumps 
disco ntinuously to the right by one unit, so that this crossing occurs at q = 1 instead of at 
q = 0. (In contrast, the right-hand part of B q increases above 2 continuously as w decreases 
below 1.) 



B. Zeros of Ph(G, q, w) in w as Functions of q 

One may also study the zeros of Ph(G,q,w) in if as a function of q. For graphs G 
containing at least one edge, Ph(G,l,w) vanishes identically We therefore take g / 1, 
although we shall consider the limit q — > 1 below. We again consider some simple examples. 
Ph(Li,q,w) = for w = 1 — q. From Eq. (13.131) . we find that Ph{L 2l q ) w) = for 
w = 1 — (q/2). From Eq. (I3.14p . it follows that Ph(L 3 ,q,w) = at w = Wl3 Z ,±(q), where 



1 

W L 3z,± = - 



5-3g± v / (g-l)(5g-9) 1 . (5.5) 



These roots are real for q < 1 and q > 9/5, and form a complex-conjugate pair for 1 < 
q < 9/5. For q — 0, these roots are 1 and 4. As q increases from to 1, u>L3 2 ,+ decreases 
monotonically from 4 to 1, while w^ z ^ first decreases, reaching a minimum of 4/5 at q = 4/5, 
and then increases to 1. Thus, at q — 1, the roots coalesce to form a double root. As q 
increases above 1, they split apart to form a complex-conjugate pair, with the magnitude of 
the imaginary part reaching a maximum at q — 7/5, for which u>L3 2j ± = (l/5)(2±v5z). As q 
increases further, these roots move back to the real axis, coming together again at w — —2/5 
as q increases through the value 9/5. As q increases further, ^l3z,- decreases monotonically, 
while wl3z,+ first increases to the value at q = 2 and then decreases. For the complete 
graph K n , Ph(K n , q, w) has a single zero in w at 

WRnz = 1 - - . (5.6) 

n 

The zeros of Ph(G,q,w) in w as a function of q are not, in general, bounded, even for 
finite values of q. This is a consequence of the fact that the coefficient of the term in 
Ph(G, q, w) of highest power of w may vanish as a function of q, in contrast to the fact that 
the coefficient of the term in Ph(G, q, w) of highest power in q is a w-independent constant 
(namely, 1) and hence never vanishes. Generically, in the absence of cancellations, a root of 
an algebraic equation diverges when the coefficient of the term of highest degree vanishes. 
This is evident in the quadratic equation aw 2 + bw + c = 0, where a, b, and c are functions 
of q with no common factors. Let us denote the set of values of q where a(q) = as {qo}- 
One of the roots of this equation, w = (2a) _1 (— b± \/b 2 — 4ac ), diverges when q approaches 
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one of the values in the set {qo}, since a — > in this limit. A similar comment applies for 
algebraic equations of higher degree. Given that we have restricted ourselves, with no loss of 
generality, to connected graphs G, these all contain at least one edge, except for the case of a 
single vertex. Hence, for these graphs with at least one edge, Ph(G, q, w) contains the factor 
(q — 1). It is thus convenient to discuss the reduced coefficients /3a,j defined in Eq. ( I2.36p . 
A zero of Ph(G, q, w) in w has a magnitude that generically diverges when the coefficient 
of the term of highest degree in w, namely the coefficient Pg,cI w (g), vanishes. This type of 
divergence can be absent if coefficient (s) Pqj with j < d w (G) also vanish sufficiently rapidly 
as q approaches the value where @G,d w {G) vanishes. 

We give some simple examples of the divergences in zeros of Ph(G, q, w) in w as a function 
of q. For the line graph L 4 , using our result in Eq. (13.151) above, we find that Ph{L^q y w) 
has zeros in w at 

w L iz,i = 2-q (5.7) 

and 

(a - 2) 2 

As q — (4/3) — ► ± , wlaz,2 =Foo. This divergence is a consequence of the fact that (i) 
the term in Ph(L±, q,w) of highest power in w, namely (q — l)(3q — 4)w 2 , has a reduced 
coefficient /3l 4 ,2 = 3g — 4 that vanishes at q — 4/3 and (ii) the terms of lower degree in w do 
not vanish at q = 4/3, as is clear from their reduced coefficients /?l 4 ,i = 2(g — 2)(2q — 3) and 
P u ,o=(q-2) 3 . 

Another example is provided by the line graph L 6 . From Eq. flB4p it follows that 
Ph(L 6 , q, w) has zeros in w at 

and 

(9-2) 

W L 6z,j = 



5 - 4q ± W8q 2 - 16q + 9 



4(,-l) • J = 2 ' 3 ' (5 ' 10 > 

where j = 2 and j '• = 3 correspond to the + and — signs, respectively. As q — (3/2) — > ± , 
w L6z,i(q) Too, and as q — 1 — > ± , %2 i2 — > =Foo while w i6 ^ 3 — > 1. The divergences 
in wlgz,i and wl§ z ,2 are consequences of the fact that (i) the term of highest power in w in 
Ph(L 6 ,q,w), namely 2(q — l) 2 (2q — 3)w 3 , has a reduced coefficient /3l 6j3 = 2(q — l)(2q — 3) 
that vanishes at q — 1 and q = 3/2 and (ii) the terms of subleading degree in w do not 
vanish at either of these values of q, as is clear from their reduced coefficients /?l 6 ,2 — 
(q - 2)(10g 2 - 28q + 19), ^ 6il = 2(3g - 4)( 9 - 2) 3 , and ^ 6j0 = (q - 2f . 
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For the I5 graph, using our result in Eq. ( 1B2j) for Ph(Y^, q, w), we find that this polynomial 
has zeros in w at 

w Y 5z,i = 2-q (5.11) 

and 



W Y5 



- 4q 2 + 13q - 11 ± (q - 1) ^8q 2 - 28q + 25 

2(2^ L - ^ 2 ' 3 ' < 5 - 12 > 

where j = 2 and j = 3 correspond to the + and — signs, respectively. As q — (3/2) — > ± , 
W5z,3 — *■ Too, while wys^ — > 1/4. These divergences can easily be understood from the 
structure of the reduced coefficients /3y 5 j for the various powers of w in Ph(Y 5 ,q,w), as 
discussed in general above. 

As a last example, for the IsoYq graph, using our calculation in Eq. flB6j) . we find that 
Ph(IsoY 6 ,q,w) has a double zero at 

wi soY 6z,i = 2-q (5.13) 

and two other zeros at 



WlsoYQ 



z,3 



- 2q 2 + 6q - 5 ± (q - l)^ 2 - 10g + 9 

- L , J = 2, 3, (5.14) 

q-2 

where j = 2,3 for the ± signs, respectively, as before. As q — 2 — > O^ 1 , w/ so y6 Z ,3 — > Too, 
while w IsoY 6,2 -> 0. 

As these examples show, the zeros of Ph(G, q, w) in u> can be unbounded as functions of 
q. The divergences in these zeros that we have found occur at values of q > 1 (including the 
integers q = 1 and q = 2). It is of interest to determine a region in the q plane for which the 
zeros of Ph(G, q, w) in w are bounded, and we are studying this problem. Our results show 
that such a region would have to exclude the real interval 1 < q < 2. 



VI. QUANTITIES DEFINED IN THE LIMIT n{G) -> 00 

A. <I> Function 

From the chromatic polynomial P(G,q) = Ph(G,q, 1), one defines a configurational de- 
generacy, which is the ground-state degeneracy, when viewing P(G, q) as the partition func- 
tion of the zero-temperature Potts antiferromagnet, 

W({G} iq )= limP(G,g) 1 /", (6.1) 
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where n = n(G) and we use the symbol {G} to denote the limit n — > oo for a given family 
of graphs (and the symbol W should not be confused with the variable w). In the present 
context, this n — > oo limit corresponds to the limit of infinite length for a strip graph of fixed 
width and some prescribed boundary conditions. The associated configurational entropy per 
vertex (ground-state entropy per site of the Potts antiferromagnet) for {G} is 

S = k B \nW. (6.2) 

(There should not be any confusion between W and w, which refer to different quantities.) 
The third law of thermodynamics states that the entropy per site S goes to zero as the 
temperature goes to zero. However, there are a number of exceptions to this law. Elementary 
lower bounds on W are W > (q— l) 1 / 2 on a bipartite graph and W > (q — 2) 1//3 on a tripartite 
graph. Hence W > 1, and S > for (i) q > 1 and (ii) q > 2 on the n — > oo limit of a 
(i) bipartite and (ii) tripartite graph, respectively. In each of these cases, the third law is 
violated. A well-known violation in nature is water ice, which exhibits ground-state entropy 
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In the present case with a nonzero external magnetic field H ^ 0, we define an analogous 
quantity 

${{G},q,w) = lim Ph{G,q,w) 1/n . (6.3) 

n(G)—» oo 

As before (cf. Eq. (1.9) of [2j| and Eq. (2.8) of 29]), one must take account of a noncom- 
mutativity of limits, namely the fact that for certain special values of q, denoted {q s }, the 
limits n — * oo and q — > q s do not commute: 

lim lim Ph(G, q, w) 1/n ^ lim lim Ph(G, q, w) 1/n . (6.4) 

n.^oo q—*q s q^qs n— >oo 

Because of these noncommutativities, the formal definition (16.31) is, in general, insufficient 
to define $ at these special points; it is necessary to specify the order of the limits that one 
uses in eq. (16.41) . This noncommutativity also affects the resultant accumulation sets B. We 
have discussed this in detail before in the case of the chromatic polynomial 



28] and zero-field 



Potts model partition function 29J. Modulo this subtlety, it follows from Eqs. (II. ip and 
(TO]) that 

9{{G},q,l) = W({G},q) (6.5) 

and 

$({G},q,0)=W({G},q-l) . (6.6) 
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B. Accumulation Locus B for Strip Graphs and $ Function 



For a n- vertex graph G in a recursive family of graphs such as strip graphs, as n — > oo, 
a subset of the zeros of Ph(G,q,w) merge to form a locus B. For fixed w, this is a locus 
B q in the g plane, while for fixed q, it is a locus £><„ in the w plane. We define a dominant 
(maximal) eigenvalue \ ma x as an eigenvalue whose magnitude |A maa; | is larger than or equal 
to the magnitudes of all other eigenvalues. From Eq. (12.301) . it follows that the zeros of 
Ph(G, q, w) can occur either as an isolated zero of a single dominant eigenvalue or where 
two dominant eigenvalues are equal in magnitude. The continuous accumulation set of the 
zeros of Ph(G, q, w) in a given variable, denoted B, is given generically by the solution set of 
the condition of equality of dominant eigenvalues. For real w, the coefficients of the terms in 
Ph(G, q, w) are also real (actually integers, although this is not used here), and consequently, 
the set of zeros of P(G, q, w) in the complex q plane is invariant under complex conjugation 
q — > q*. A fortiori, the locus B q is also invariant under this complex conjugation. By the 
same logic, for real q, the set of zeros of Ph(G, q, w) in the complex w plane is invariant 
under the complex conjugation w — > w*, and so is their accumulation set B w . Because the 
A's are the same for lattice strip graphs with cyclic and Mobius boundary conditions, it 
follows that in the L x — > oo limit, the loci B are also the same. With regard to the zeros of 
Ph(G, q, w) in w for fixed q, we note that these are to be contrasted with the zeros of the 
Potts model partition function Z(G, q, v, w) in w for fixed q and v ^ —1, which have studied 
previously in many works. 

For the strip graphs of width L y considered here, 

mG},q,w) = (\ max y/ L y . (6.7) 

As one moves across a locus B, there is thus, generically, a switching of dominant eigenvalues 
and an associated non-analyticity in $. From Eq. (11.21) we have 

$({G},g,0) = $({G},g-l,l) . (6.8) 

VII. $ FUNCTION AND ACCUMULATION LOCUS B FOR LINE GRAPHS 

A. B q 

Only Ai i0 j, j = 1,2 contribute to Ph(L n ,q,w). For w = 1 we encounter the noncom- 
mutativity of Eq. (16.41) . If we first set w — 1 and then vary n, we can use the fact that 
Ph(L n , q, 1) = P(L n , q) = q(q — l) n_1 , so that aside from the single zero at q = 0, the zeros 



35 



accumulate at q = 1 and B q degenerates to this single point. If we choose the other order of 
limits, first taking n — > oo and then w — > 1, then the locus B q is the solution to the equation 
| A i ,o, 1 1 = 1-^1,0,2 1- This equation is always satisfied if q = 2, so this point is on B q , which forms 
a complex-conjugate arc with endpoints where A\ = 0, where A\ was given in Eq. (I4.23p . 
As w — > 1, these endpoints come together at q — 0. 

For w = 0, Ph(L n , q, 0) = P(L n , q — 1) = (g — l)(g — 2) n_1 , so the locus B q degenerates 
to the single point at q — 2. We proceed to consider w ^ 0, 1. Here, the equation 
| A 1 5 o, i I = I ^1,0,2 1 determines the locus B q . The DFCP interval < w < 1 is of particular 
interest, since for this interval the weighted chromatic polynomial Ph(G,q,w) interpolates 
between two chromatic polynomials; Ph(G,q, 1) = P(G,q) and Ph(G,q,0) = P(G,q — 1). 
We thus study B q for this interval first. In this DFCP interval, B forms a (self-conjugate) 
arc passing through q = 2 and ending at the points where A\ = 0, namely q e j, j = 1, 2, 



l-w± ^w(w - 1) 1 , j = l,2, (7.i; 



where j = 1,2 correspond to the ± signs, respectively. Here the square root in Eq. ( 17. ip is 
pure imaginary, so q e ^\ — q* 2 form a complex-conjugate pair. As w decreases below 1, these 
arc endpoints move away from the real axis near q = 0. They reach their maximal distance 
from the real axis at w = 1/2, where q e j = 1 ± i, and as w decreases from 1/2 to 0, these 
arc endpoints come back toward this axis, finally reaching it at q — 2. 

In the FCP interval w > 1, the locus B q is a line segment whose right and left ends occur 
at q e ^i and g ej 2, respectively. As w increases above w — 1, the line segment extends outward 
from the point q = 0. Some illustrative sets of values for these endpoints are (i) for w = 1.2, 
q eA ~ 0.580 and g e , 2 ~ -1.380, (ii) for w = 2, q e>l ~ 0.828 and g e , 2 ~ -4.828, and (iii) for 
w = 10, q e> \ = 0.974 and q £)2 = —36.974. As w — > 00, the the right end of this line segment 
occurs asymptotically at 

qe,i = l-^--^-0(\) , (7.2) 

4w 8«r V-ur/ 

while the left end occurs at approximately 

q e . 2 = -Aw + 3 + ^- + (\) . (7.3) 
Aw \w z / 

As discussed above, the ranges of w for these weighted coloring problems are < w < 1 
and w > 1, respectively. We may also consider an extension of this range of values of 
w in which w becomes negative, although negative values of w do not correspond to a 
weighted graph coloring problem. As w decreases from through negative values, B q forms 
a line segment that extends outward from the point q = 2, with left and right ends at 
g ei i and q e ^ respectively. As w — > —00, the left end approaches q = 1 from above, as 
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1-1 



g 6j i ~ l + (A\w\)~ 1 + 0(w~ 2 ), while the right end goes to infinity, as q e ^ — 4| w| + 3 — (4|w| 
This locus B q does not separate the complex q into any separate regions. The dominant A 
is Ai.,0,1 an d, denoting the formal limit of this family of line graphs as lim^oc L n = {L}, the 
resultant $ function is 



$({L},q,w) = Ai.o.1 



1 



-2+ v / (g-2) 2 + 4(g-l) 



w 



(7.4) 



Here it is understood that one takes account of the branch cut associated with the branch 



point singularity in the square root, so that at large negative q, | A i ,o, 1 1 
this has the Taylor series expansion 



-q. For w ~ 1, 



®({L},q,w) = q- 1 



-!)(«,-!)] [(q - l)(w - l)f 



+ 



2[(g-l)( W -l)P 



as w 



For w ~ 0, the Taylor series expansion for q, w) is 



$({L},g,u7) = (g-2) 



1 + z - z 2 + 2z 3 - 0(z 4 



as w 



where here we use the compact notation 

(g - l)w 

2 (g-2) 2 ' 

oo, w) behaves asymptotically as 

q-2 



(7.5) 



(7.6) 



(7.7) 



As lid 



${{L},q,w)~y/(q-l) 



w 



1 + 



2y/(q-l)w ^wJ 



as 



oo 



(7.1 



As g 



oo, $({L},g, w) behaves asymptotically as 

w(w — 1) 2w(w 



$({L},q,w) ~q-2 + w 



+ 0[ — ) as q — > oo 



(7.9) 



In Fig. [2] we show plots of $({L}, q, w) as a function of w for some representative values 
of q. As is evident from Fig. [2J for g > 1, $({L}, g, w) is a monotonically increasing function 
of w in the DFCP range < w < 1. This reflects the fact that it is easier to carry out a 
proper g-coloring of the line graph as the weighting factor w increases from to 1. It is also 
easy to understand why $({L},g, w) is, for g > 1, a monotonically increasing function of 
w in the FCP interval of real w > 1, with the leading asymptotic form given in Eq. (17. 8p . 
$({L}, g, w) ~ y(q — T)w. To see this, let us start with finite n and maximize the number 
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FIG. 2: Plot of <J>({L},g, w) as a function of w for the following values of q, from bottom to top: 
q = (a) 2 (b) 2.5 (c) 3 (d) 3.5 (e) 4. 

of vertices assigned the color 1, in order to maximize the power of w in $(L n ,g, w). Since 
the graph L n is bipartite, we can start at, say, the left end of the line graph and assign the 
corresponding vertex, denoted vertex i = 1 with this color 1. Then the next vertex to the 
right, i = 2, cannot be assigned this color, but there are q — 1 possibilities for its color. 
The vertices i = 3, 5, and so forth for all odd-numbered vertices, are similarly assigned the 
color 1. Each of the even-numbered vertices can independently be assigned any of q — 1 
colors. Taking into account all of these possible color assignments (or equivalently, spin 
configurations, in the statistical mechanics context), it follows that, if n is even, then the 
dominant term in &(L n ,q,w) as w — > oo for q > 1 is Ph(L n , q,w) = [(q — l)w] n ^ 2 and if 
n is odd, then this dominant term is Ph(L n , q,w) — (q — l)( n_1 )/ 2 u;( n+1 )/ 2 . In either case, 
for q > 1, as w — > oo, if one takes n — > oo and calculates $({L},g,w), one obtains, as the 
leading asymptotic expression, $({L},g, w) ~ \/{q — l)w. 

Although negative w is not associated with any coloring problem, we also show in Fig. [2] 
the extensions of the curves into the negative-u> region. The function A\ in the square root 
of $({L},g, w) becomes negative, and hence w) becomes complex, for w < w z (q), 

where 

(a - 2) 2 

•"■<«> = -iH) ' (7 ' 10) 

For example, w z {2) = 0, w z (3) = —1/8, and u> 2 (4) = —1/3. We only plot each curve for 
w > w z (q). The values at w = and w = 1 are $({L}, 0, w) = q — 2 and $({L}, l,w) — q — 1. 
In Figs. [3]andH]we show q, w) as a function of q for the range relevant for coloring, 
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V ' ' 1 .5 ' ' ' 2 ' ' 2.5 ' ' ' 3 ' ' ' 3.5 ' ' ' 4 
1 



FIG. 3: Plot of <£({L}, q, uu) as a function of q > 1 for the following values of w: w = (a) 0.1 (b) 
0.2 (c) 0.5 (d) 0.8 (e) 0.9. The plot also shows Ph{{L},q,0) = q - 2 and Ph({L},q, 1) = q- 1. 
The curves for the w values (a)-(e) are arranged from bottom to top, to the right of q = 1, between 
these lines. 




FIG. 4: Plot of <1?({L}, g, as a function of g > 1 for the following values of w in the FCP interval: 
w = (a) 1.1 (b) 1.5 (c) 2 (d) 3. The curves for the w values (a)-(d) are arranged from bottom to 
top. 
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namely q > 1, and a set of w values in the DFCP interval < w < 1 and the FCP interval 
w > 1, respectively. One sees that for fixed g > 1, $({L}, g, w) is a monotonically increasing 
function of w for w > 0. Given its definition as the 1/n'thd power of the weighted chromatic 
polynomial Ph(G, q, w) as n — > oo and the fact that Ph(G, q, w) provides a measure of 
the ease of carrying out a weighted proper g-coloring of the graph G, it follows that where 

w) is meaningful for such colorings, it is non-negative. Hence the contination of 
the line that passes through q = 2 to negative values of $ is not relevant for graph coloring. 

These calculations with line graphs show a number of general features of the weighted 
chromatic polynomial Ph(G,q,w) and the associated limiting function <&(G,q,w). They 
suggest the following conjectured generalization for weighted colorings of families of strip 
graphs G of regular lattices A: in the n — > oo limit of G, denoted {G}. Let the chromatic 
number of the lattice A be indicated as For a bipartite lattice such as a line graph 

L n , a circuit graph C n with even n, or a square, cubic, or body-centered cubic lattice, 
x(A-bip.) = 2. For a triangular lattice, x(^ r — 3, etc. Assume q > Then (i) for 

fixed w > 0, &({G},q,w) is a monotonically increasing function of q and (ii) for fixed q, 

w) is a monotonically increasing function of w for w > 0. The generalization (i) 
is understandable since q, w) is a measure of the number of weighted proper g-coloring 
of G, and this should increase if there are more colors, i.e., if q increases. For the DFCP 
interval < w < 1, the generalization (ii) follows because increasing w in this interval 
removes the penalty factor for coloring a vertex with one of the colors and hence clearly 
makes it easier to perform a proper g-coloring of the vertices of G. In the FCP interval 
w > 1, the generalization (ii) is rendered plausible because for sufficiently large q, one can 
analyze the dominant terms contributing to <&(G,q,w), and these arise from maximizing 
the number of vertices that can be assigned the color 1, subject to the proper g-coloring 
condition, and then enumerating the possible color assignments for the other vertices. We 
remark that for these monotonicity properties, it is important that q > x{G)- As an example 
of how the behavior differs when q < x(G), consider the weighted chromatic polynomial of 
the complete graph, K n , given in Eq. (13.241) . Recall that x{Kn) — n - Let us take n = 4 for 
definiteness, so that Ph(K 4 , q, w) = (q — l)(q — 2) (5 — 3)[q+4(w — 1)]. This is a monotonically 
increasing function of q if q > 3, but not for smaller values of q. Moreover, say we keep w 
arbitrary but choose the illustrative value q = 5/2, whence Ph(K±, 5/2, to) = (3/16) (3 — 8w). 
This is not an increasing function of w in the DFCP or FCP intervals. 
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One can also study B w as a function of q. We find that B w is the semi-infinite real line 
segment 

B w : w<w z (q) for {L} , (7.11) 

where w z (q) was given in Eq. (17.101) . For the range of q relevant for weighted graph coloring, 
namely q > 1, w z (q) < 0. 



VIII. $ FUNCTION AND ACCUMULATION LOCUS B q FOR CIRCUIT 

GRAPHS 

Here we discuss the accumulation set B of the zeros of Ph(G,q,w) as n — > oo for the 
family of circuit graphs, C n (this limit is denoted {C}). The results depend on w, so we 
discuss various intervals of w in turn. 



w = 1 



We begin by briefly reviewing the results for the case unweighted case w = 1, i.e., for 
the chromatic polynomial P(C n , q), given in Eq. (14.301) . The resultant locus B q in the limit 
n —>■ oo is the unit circle 



B q : \q — 1| = 1 for w = 1 . 
This crosses the real axis at the points 

q cr ,i = , q cr ,2 = 2 for w = 1 



U) 



(8.2) 



Following our earlier notation 28|, |30j, |31( , we denote the maximal point at which B q intersects 
the real q axis for the n — > oo limit of a given family of graphs, {G}, as q c ({G}). Thus, here, 
q c = 2. 



Indeed, all of the complex zeros lie exactly on this unit circle 32]. The boundary B q 
separates the q plane into two regions, in which W(q) has different analytic forms. Outside 
of the circle \q — 1| = 1, the dominant A is A^o.i — q — 1, while inside of this circle, the 
dominant A is Ai 5 i = —1, so 

W(q) = q - 1 for \q - 1| > 1 



\W(q)\ = 1 for \q- 1| < 1 . 



(8.3) 
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B. < w < 1 



As w decreases below 1 in the interval < w < 1, the boundary B q continues to be a 
simple closed curve separating the q plane into two regions. However, there is a discontinuous 
change in the form of this boundary. On the left, the point at which the locus B q crosses the 
real q axis jumps from q = for w = 1 to 

g cr ,i = 1 for < w < 1 . (8.4) 

Associated with this, the left part of the boundary B q changes discontinuously from a section 
of a circle to an involuted cusp with its tip at q — 1. (There may also be discrete zero(s) to 
the left of q = 1.) Thus, as one moves along the curve forming B q upward from the point 
q = 1 where it intersects the real axis in the cusp, this curve moves to the upper left, finally 
curving around to go upward, and then over to the right. The behavior of the boundary B q 
on the right side is continuous as w deviates from 1; this part of B q crosses the real axis at 

q C r,2 = q c = ^— for {G} = {C} . (8.5) 
w + 1 

This point q c increases monotonically from q c = 2 for w = 1 to q c = 3 as w decreases from 
1 to 0. We denote as Ri the region that includes the real interval q > q c and the part of 
the complex q plane analytically connected to it, which is the region outside of the closed 
curve formed by B q , For w only slightly less than unity, B q has the form of a lima bean, 
with its concave part facing left and its convex part facing right. As w decreases through 
this interval < w < 1, the bulbous parts of B q on the upper left and lower left disappear, 
and eventually, as w decreases toward 0, the locus B q becomes the circular locus (18.1 p with 
q replaced by q — 1 (in accordance with Eq. (II. 2p ). i.e., unit circle whose center is shifted 
horizontally by one unit to the right in the q plane: 

B q : \q - 2| = 1 for w = . (8.6) 

In region R±, 

<S>({C},q,w) = mL},q,w) for q G R 1 , (8.7) 

where $({L},g, w) was given above in Eq. (17.41) . In region i?2 forming the interior of the 
closed curve B q , 

m{C},q,w)\ = 1 for qeR 2 - (8.8) 

Thus, a plot of <3>({C},g, w) as a function of q for w G (0,1) is similar to the plot of 
$({L},g, w) given in Fig. [3] but with the difference that the curve extends only down to 
the value q = q c given in Eq. (18. 5p . where $({L},g c ,w) = $({C},g c ,w) = 1, and for 
1 < q < q c , &{{C}, q, w) has unit magnitude (in region R 2 , while <3>({L},g, w) continues 
downward, reaching zero as q — > 1. 



42 



1. w > 1 



For the range w > 1, the locus i3 g contains a line segment on the real axis, whose left 
end occurs at q e ^ (cf. Eq. (17.10 ). This is a solution of the condition that A\ = 0, where A 1 
is the function in the square root in Eq. ( 14.231) . Along this line segment, A\ < 0, so that 
the square root in Eq. (14.221) is pure imaginary; hence, the Xipj with j — 1, 2 are complex 
conjugates of each other. They are also larger than |Ai,i| = 1 and hence this line segment is 
on B q . As one moves to the right along the real axis on this line segment, when one comes 
to the intermediate point q int , given by 

Qint = 1 - - , (8.9) 
w 

the equal magnitudes of Xiflj, j = 1,2 decrease through the value 1, so this is the point at 
which this line segment on B\ terminates. For the present range w > 1, the point q = is 
always on this line segment, since 



Ai,o,i(<? = 0) = -l±*vV^I, j = l,2, (8.10) 

and these are dominant over A^i = —1, since IA^qjI = 1 + \w — 1| for j = 1,2. In the real 
interval 

q in t<q<l, (8.11) 

and in the region of the complex q plane analytically connected with it, Ai ; i = — 1 is domi- 
nant. We denote this region as i? 3 , and 

\$({C},q,w)\ = l for qeR 3 . (8.12) 

Thus, although the square root in the Xi t o,j is imaginary for the full interval q e ^ < q < q e ,i, 
this does not affect B q for q > q int . 

At q — 1, Ai i0 ,2 becomes degenerate with A^i, so B q crosses the real q axis at this point. 
The eigenvalue \i,q,2 is dominant in the interval 

l<g<2 (8.13) 

and the region of the q plane analytically connected with it. We denote this region as R2, 
and obtain 

$({C}, q,w) = - q - 2 - y/fa - 2) 2 + 4{q - l)w for q G R 2 . (8.14) 

At the point 

q = q c = 2 (8.15) 



43 



-4 




FIG. 5: Plot of (a) | Ai o,l I j 3 ) I ^1,0,2 1) an d (c) \Xi i\ = 1 for Ph{C n ,q,w) with the illustrative value 
w = 3/2. The values of these magnitudes and the corresponding order of the curves from top to 
bottom are (i) a > b > c for q > 2; (ii) b > a > c for q E (1, 2), (iii) c > b > a for qi nt < q < 1, (iv) 
a = b > c for q ei 2 < q < qint, (v) a > b > c for q < g e ,2- For if = 3/2, g in j = 1/3 and q e> 2 — —2.732. 
See text for further discussion. 

Ai,o,i = — ^1,0,2 = and both are dominant over A^i = —1, since w > 1. Hence, £> 9 
crosses the real q axis again at this point, and this is the maximal real value of q where B q 
crosses the real axis. We denote the real interval q > 2 and the region of the complex q plane 
analytically connected with it as region R^. In this region, since Ai o,i is dominant, 

q,w) = ^q-2 + y/(q- 2) 2 + A{q - l)w ] for q G R x . (8.16) 

Thus, for w > 1, the locus B q for the n — > 00 limit of Ph(C n , q, w) separates the g plane into 
three regions. 

In Fig. Owe plot |Ai,o,i|, |Ai,o,2|, and [ Ai^i | = 1 as functions of q for the illustrative value 
w = 3/2. For this value, q int = 1/3 and g 6i2 = —(1 + ) ~ —2.732. As before, the A of 
dominant magnitude determines <&({C},q,w). For q > 2, Ai o,i is dominant; for 1 < q < 2, 
Ai,o,2 is dominant; for 1/3 < q < 1, Ai,i is dominant, and for q < g e 2 , Ai^i is again dominant 
(taking into account the branch cut in the definition of the sign of the square root). On the 
interval q e ^ < q < qint, |Ai,o,i| and ^1,0,2 1 are equal in magnitude and dominant. This plot 
shows the intersections of the curves (or line) where there is a degeneracy of dominant A's 
at q — 2, q — 1, and along the interval q e $ < qint- 



44 



C. w < 



As w decreases below through the range — 1/3 < u> < 0, B q forms a closed curve that 
crosses the real axis at the fixed point q = q CTy \ = 1 and at the w-dependent point q = q cr ^ 
which increases as w becomes more negative. The shape of B q changes as w becomes more 
negative in this interval, becoming a teardrop with its broadly rounded end on the left and 
its sharper end on the right. As w decreases through the value w = —1/3, a line segment 
appears, and for w < —1/3, the rightmost part of the locus B is comprised by this line 
segment, which extends from q = q int in Eq. (18.91) to q e> i in Eq. (17. ip . Thus, for w < —1/3, 
the locus B consists of a teardrop-shaped curve crossing the real axis on the left at q — 1 
and on the right at q — qi nt (with the bulbous part of the teardrop on the left and the sharp 
part on the right), together with a line segment extending over the interval q^ < q < q Bt \. 
Mw = -1/3, 

(Ai, .iU=-i/3 = ~ [(q ~ 2) ± v/(g-4)(g-(4/3)) ] . (8.17) 

For this value w = —1/3, and q = 4, the pair Xi,o,j, j = 1,2, are equal to each other 
and are also equal to the magnitude [ Ax,i | = 1- As w — > —oo, q int — > 1 + , so the teardrop 
curve contracts to a point at q = 1 while the right-hand endpoint of the line segment at g e l 
approaches infinity like q e ^ ~ 4|u>|. 

One can also calculate the loci B q for the n — > oo limits of other families of graphs. 
However, our discussion for the n —>■ oo limits of line graphs L n and circuit graphs C n 
already exhibit a number of salient features of these loci. 



IX. LOCUS B w FOR CIRCUIT GRAPHS 

Here we discuss the locus B w as a function of q for the n — > oo limit of the circuit graph 
C n . First, we note that if q = 1 or q = 2, then we encounter the noncommutativity (16. 4p . 
For q — 1, this is evident from Eq. (12.181) . according to which if we set q — 1 first and then 
take n —>■ oo, the problem is trivial, since Ph(C n , l,w) vanishes identically. For q = 2, the 
noncommutativity is evident from our general result in Eq. (I4.32p . because the coefficient of 
the (Xi,i) n term in Ph(C n , q, w) vanishes if q — 2. If we first take n —>■ oo and then set q = 2, 
we have Ai^j = ±y/w, so that the locus B w is the union of the semi-infinite real line segments 
w > 1 and w < — 1. If, on the other hand, we first set q = 2, and then vary n, our result 
(I4.32p shows that the limit of Ph(C n ,2,w) as n — > oo does not exist, since Ph(C n ,2,w) is 
alternatively zero for odd n and 2w n ^ 2 for even n. If we restrict to odd n, then the problem 
of the zeros of Ph(C n , 2, w) is trivial since the function itself vanishes identically, while if we 
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restrict to even n, then the locus B w degenerates to a point at w = 0, correswponding to the 
zero at this point with multiplicity n/2. 

From Eq. (12.221) we know that for q = 0, Ph(G,0,w) contains a factor of (w — 1) for 
an arbitrary graph G, and this is true, in particular, for G = C n . The other zeros occur at 
real values w > 1. For this value q — 0, it follows that Xi t oj = — 1 ± y/\ — w, and these are 
dominant A's if w ^ 1, so B w is the semi-infinite real line segment w > 1. For g ^ 0, 1, 2 
we typically find that the locus B w may consist of the union of a (self-conjugate) loop and a 
line segment. Details depend on the specific value of q. 



X. B FOR WHEEL GRAPHS 

We have also calculated the locus B for wheel graphs. We denote the n — > oo limit of the 
graph Wh n as {PF/t}. For w — 1, B q is the unit circle |? — 2| = 1, which crosses the real axis 
at q = 1 and g c = 3 and separates the q plane into two regions. In the region with \q — 2\ > 1, 
i.e., the region exterior to this circle, §({Wh}, q, 1) = W({W7i}, q) = q — 2. In the for which 
\q — 2| < 1, i.e., the region interior to the circle, q, 1)| = q) \ = 1. 

As u> decreases from 1 in the DFCP interval, the boundary B q continues to form a closed 
curve separating the q plane into two regions, as it did for w — 1, but there is a discontinuous 
jump in the crossing point on the left, from q — 1 to q = 2. This is similar to what we found 
for the n — > oo of the circuit graph, where the jump in the crossing point on the left was 
from q = to q = 1. The crossing point on the right is 

g c = ^^ for {G} = {C} (10.1) 

Region i?i includes the real interval q > q c and the portion of the complex q plane analytically 
connected with this interval, and thus lying outside of the boundary B q . Region R 2 occupies 
the portion of the q plane inside of the boundary B q . The dominant A in region Ri is Xwh,+, 
while the dominant A in R 2 is equal to —1. The point q c occurs where these are degenerate in 
magnitude. Given this and the relation (13.331) . it follows that q c for {G} = {Wh} is related 
to q c for {G} = {C}, given in Eq. (18.51) . by replacing q by q — 1. That is, if one replaces q c 
on the left-hand side of Eq. (18. 5p by q c — 1 and solves for the new q c , one obtains Eq. (110. ip . 
This is in accord with the fact that the proper g-coloring of the wheel graph with q colors is 
closely related to the proper coloring of the circuit graph with q — 1 colors. As w decreases 
from 1 to in the DFCP interval the q c in Eq. (UO.ip increases continuously from 3 to 4. 
One can also analyze other ranges of w and the locus B w in a similar manner. 
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XI. SOME OBSERVATIONS AND CONJECTURES 



A. Sign Alternation of Successive Terms in Ph(G, q, w) 
One can write the chromatic polynomial P(G, q) of a graph as 

n-l 
3=0 

where, without loss of generality, we take G to be connected. The signs of the coefficients 
ac,n-j alternate: 

sgn(« G , n _,) = , < j < n - 1 . (11.2) 

This is proved by iterated application of the deletion-contraction theorem. Since the weighted 
chromatic polynomial Ph(G, q, w) does not, in general, obey a deletion-contraction theorem, 
except for the values w = 1 and w = for which it reduces to a chromatic polynomial 
(see Eqs. (11.11) and (11. 2p ). one does not expect the coefficients ac,n-j(w) in Ph(G,q,w) to 
have this sign-alternation property, and they do not. However, from our analysis of weighted 
coloring polynomials for several families of graphs, we have noticed that for a restricted range 
of w, namely < w < 1, this sign alternation again holds, namely sgn(aG,n-j(w ; )) = ( — l) 5 
for < j < n — 1. For j = n, namely for the q° term in Ph(G,q,w), the sign alternation 
also holds for < w < 1; here the coefficient ac,o contains the factor (w — 1) and hence 
vanishes at w — 1. We also find that this sign alternation property holds, as far as we have 
checked it, for real negative w. It is of interest to investigate whether this sign alternation 
property for w < 1 holds on other families of graphs. We are currently continuing with this 
investigation. 

Related to this, it is of interest to study where the coefficients otG,n-j(w) vanish in the 
complex w plane. From our calculation of weighted chromatic polynomials for line and 
circuit graphs L n and C n , we have observed that for the graphs we have considered, the 
coefficients «L„,n-i and ac n , n ~j for 1 < j < n — 1 have zeros in the real interval w > 1, 
while for the coefficients ai n ,o and ac„,o have, in addition to the always-present zero at 
w = 1 (recall Eqs. f!2.22j) and (12.331) ). the other zeros, if any, again occur in the real 
interval w > 1. In contrast, we find that for other graphs, the coefficients otG,n-j may 
have complex-conjugate pairs of zeros. For example, in Ph(S^ q, w), the coefficient of the 
q term, a5 4i i = w 3 — 9w 2 + 27w — 20, has zeros at w ~ 1.087 and w ~ 3.9565 ± 1.6566z, 
and the coefficient of the q° term, as 4 ,o — —( w — l)(w 2 — 5w + 8), has zeros at w = 1 and 
w = (l/2)(5±\/7 i). Similarly, coefficients of star graphs S n for larger n include cases having 
complex-conjugate pairs of zeros in the w plane. 
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B. Generalized Unimodal Conjecture 

From his study of chromatic polynomials, R. Read observed that the magnitudes of the 
coefficients of successive powers of q n ~\ < j < n — k(G) in a chromatic polynomial satisfy 
a unimodal property |8]. That is, the magnitudes of these coefficients get successively larger 
and larger, and then smaller and smaller, as j increases from to n — k(G). There is thus a 
unique maximal-magnitude coefficient, or two successive coefficients whose magnitudes are 
equal. From our calculations of weighted chromatic polynomials for a number of families of 
graphs, we have observed that in the interval < w < 1 this property continues to hold. 
We therefore state the following conjecture: Conject. Let Ph(G,q,w) be written as in Eq. 
f)2.30p . Then for real w in the interval < w < 1, the quantities (-l) J «G,n-j(w), < j < n, 
are positive and satisfy the unimodal property, i.e., (— l)i a.G, n -j{ w ) get progressively larger 
and larger, and a maximal value is reached for a given j, or for two successive j values, and 
then the quantities (— l)%G,n-j(^) get pogressively smaller, as j increases from to n. 



XII. SOME GENERALIZATIONS 



Although we have focused in this paper on the proper g-coloring of vertices with one 
color given a disfavored or favored weighting, we discuss some generalized weighted coloring 
problems in this section. We first present an extension of Eq. (12. 4p to the most general case 



of different fields H p , p 



q, and hence different weighting factors, for each of the q 



different colors. The generalization to multiple fields corresponding to different spin values in 
the Potts model was noted, e.g., in 
for this case is 



331 ] and 34j and more recently in 25|. The Hamiltonian 



H 



EKE- 



■a e ,p 



Let us define 



h p = (3H p 



w, 



for 1 < p < q 



;i2.i) 



;i2.2) 



and denote the set of w p , p = 1, ...,q as {w}. The partition function is a function of q, v, 
and {w}, and hence we write it as Z(G,q,v, {w}). For v = —1, the resultant generalized 
weighted chromatic polynomial is Ph(G, q, {w}) = Z{G, q, —1, {w}). 

We have derived the following generalization of the Wu formula for this case (where, as 
before, G' is a spanning subgraph of G): 

fe(G') 



z( G , q ,v,{w})=j2 ve{G,) n (£< (g:) ) 



;i2.3) 



G'CG 



p=l 
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This is proved as follows. The spins in each component G\ of G' are connected by edges, so 
they all have the same value, and there are q possibilities for this value. For a given spanning 
subgraph G', the weighting factor is the product Yli=i ^ ( Y^p=i w p^ G ^ ■ This subgraph thus 

contributes a term t> e ( G ') Ylt^i ^ { Sp=i Wp^^ 1 to Z. Summing over all spanning subgraphs 
G' then yields the result il2l . □ 

The resultant spanning graph formula for the generalized weighted chromatic polynomial 
Ph(G, q, {w}) is obtained by evaluating Eq. (j!2.3p at v — — 1: 

p h (G, g , {w}) = y: (-i) e(G,) n ( e < Gd ) • ( i2 - 4 ) 

G'CG i=l p=l 

Note that some w p 's may disfavor certain color (s), i.e., < w p < 1, while others may favor 
other color (s), uy > 1. Note also that in the general situation with different H p , p = 1, q, 
the dependence of Ph(G,q, {w}) on g appears via the w p , p = l,...,q rather than via a 
polynomial dependence on the variable q. 

Let us illustrate this generalization for the case where a set of s colors is subject to a 
given (disfavored or favored) weighting, i.e. (where without loss of generality, we label these 
s colors as 1, 2, s) 

( H ^ for 1 < p < s 
H p = { (12.5) 
I for s + 1 < p < q 

so that 

{w ^ 1 for 1 < p < s 
(12.6) 
1 for s + 1 < p < q 

Then, with Z(G,q,v, {w}) written compactly as Z(G,q, s,v,w) in an evident notation, Eq. 
( 112. 3p takes the form 

fc(G0 

Z(G,q,s,v,w) = ^ w e(G " ) J"| ( q - s + sw n (12.7) 

G'CG i=l 

(To avoid awkward notation, we use the same symbol Z for the Potts model partition 
function with the various sets of arguments, Z(G,q,v), Z(G,q,v,w), Z(G,q,v,{w}), and 
Z(G,q, s,v,w).) Writing the weighted chromatic polynomial Ph(G,q,v,{w}) in the same 
notation as Ph(G, q, s, w), we have 

Ph(G, q, s, w) = Z{G, q, -1, s, to) . (12.8) 

With Yan Xu at Stony Brook, a study of the properties of the generalized weighted chromatic 
polynomial Ph(G,q,s,w) has been carried out, and the results will be reported elsewhere. 
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We note here that since s only appears in Eq. ( 112. 7ft in the combination 

fc(G') fc(G') n(G' t )-l 

n ^ n(G °) = n ( g + s ( w ~ ^ £ > (12.9) 

t=l i=l r=0 

it follows that Z(G, q, s, w) and Ph(G, q, s, w) can equivalently be written as polynomials in 
the variables q, v, 

t = s(w-l), (12.10) 

and w. We mention the following general relations involving Z(G, q, s, v, w) and 
Ph(G, q, s, w) that hold for t = : 

Z(G, q, 0, v, w) = Z(G, q, s, v, 1) = Z(G, q, v) , (12.11) 

Ph{G, q, 0, w) = Ph{G, q, s, 1) = P(G, q) , (12.12) 
and the general relations that hold for w = 0: 

Z(G,q,s,v,0) = Z(G,q-s,v) , (12.13) 

and 

Ph(G, q, 8, 0) = P(G, q - s) . (12.14) 

From Eqs. (112.111) and (112.121) . it is clear that for s — 0, the Z and Ph polynomials reduce 
to their zero-field forms. A similar reduction to a factor times the zero-field forms occurs if 
s = q: 

Z(G, q, q, v, w) = w n ( G) Z(G,q,v) (12.15) 

and 

Ph{G, q, q, w) = w n(G) P{G, q) . (12.16) 

Hence, we are primarily interested in the values s = l,..,q — 1. Assuming that G contains 
at least one edge, then, if q — 1, it is impossible to satisfy the proper g-coloring constraint, 
so Ph(G, 1, s, w) = 0. This vanishing does not, in general, result as a consequence of an 
explicit (q — 1) factor, unless s = or s — 1. Instead, when one sets q = 1 in Ph(G, q, s, w), 
one obtains a polynomial with a factor of s(s — l)(w — l) 2 . Since s is a non-negative integer 
bounded above by q, the condition that q = 1 implies that s is either or 1, and hence this 
factor must vanish, yielding the necessary result that Ph(G, 1, s, w) = 0. 

A different type of generalization is to have the spin-spin couplings depend on the edges of 
the graph G, so they would be of the form J^- = J e , where i and j denote adjacent vertices 
of G connected by the edge e. The study of spin models with spin-spin couplings that 
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are different for different lattice directions goes back to the early decades of the twentieth 
century, reflecting the fact that there are often anisotropies in real magnetic substances. In 
the 1960's and 1970's, anisotropic spin-spin couplings were studied to investigate how, for 
various ferromagnetic and antiferromagnetic combinations on different lattices, they could 
affect critical behavior 35|. In the 1970's and later the further generalization to spin-spin 

i disordered materials 



361 ] : and discussions of 



couplings that depend on each edge was studied in connection wit 
and the question of how such disorder changed the critical behavior 
edge-dependent Jy continue [25j. The Hamiltonian for the general Potts model for this case, 
including the full set of q different external fields, is 



H 



(ij) 



p=l 



Let us define 



(3 Jij, 



;i2.i7) 



'12.18) 



and denote the set of all as {v}. The partition function is then Z(G,q, {v}, {w}). We 
give the following general formula for this partition function, where again G' = (V, E') is a 
spanning subgraph of G: 



.KG') 



z(G, q ,{v},{w})= j2 [rH [n (e 

G'CG 



w 



12.19) 



e£E> 



p=l 



For the weighted proper g-coloring problem, v e = — 1 V e G E, so this generalization reduces 
to Eq. (I12.4p . Having derived and presented spanning graph formulas for the Potts model 
partition function and weighted chromatic polynomial for these generalized cases, we focus 
henceforth on the simple case where only one color is subject to the (disfavored or favored) 
weighting. Note, as before, that in the general situation, the dependence of Z(G, q, {v}, {w}) 
on q appears via the w p , p — 1, q rather than via a polynomial dependence on the variable 
q. We illustrate these generalizations with the circuit graph G = C 3 . Let us define r\ r = 
Y^l=i w p- Then from Eq. (112.191) we have 

Z(C 3 ,q,{v},{w}) = T)l+(vi2+V23+V 3 l)T]2'ni+ (^12^23+^23^31 +^31^12)+^ 12^23^31 V2, (12.20) 

and, setting v e = — 1 for all of the edges in C 3 , we obtain Ph(C 3 , q, {w}) — rjf — 3772^1 + 2r] 3 . 

Yet another generalization is to make the sets of colors that one chooses from to assign 
to each vertex depend on the vertex. With the weighting, this defines a new weighted list- 
coloring problem. A practical realization of this problem is the allocation of frequencies 
to radio broadcasting or wireless mobile communication transmitters where each individual 
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transmitter has its own set of available frequencies, no adjacent transmitters should use the 
same frequency, and there are various disfavored and/or favored frequencies. For a graph 
G = (V,E), we denote the list of available colors for a given vertex as {q}, where % = 
1, ...,n(G), and we denote the set of all color lists with the symbol {{c}} = {{ci}, {c n }}. 
We define the associated partition function as Z(G, {{c}}, {v }, {w}) = ex P(~ 
with 

q 

k = - E - E [ h p E <w] > ( 12 - 21 ) 

(ij) p=i e 

where Oi takes on values in the list {q}. As before, one may consider special cases of this 
weighted list coloring problem in which, e.g., the are constants, independent of the edge 
joining the vertices % and j, and/or where the H p are of the simple form (112.51) . etc. As a 
simple example, we again take the circuit graph G = C3 and choose the available color lists for 
each vertex as {ci} = (1, 2), {c 2 } = (2, 3), {c 3 } = (1, 3). The generalized weighted chromatic 
polynomial for this weighted list coloring problem would then be 2W1W2W3 corresponding to 
the color assignments (1, 2, 3) and (2, 3, 1) to vertices i — 1, 2, 3. In passing, we recall the 
case where there are no external fields or corresponding weightings, i.e., H p = 0, so w p = 1 for 
all p — 1, q. This is the usual unweighted list coloring problem, as reviewed, e.g., in Ref. 



371 ] . For example, for the case G = C3 with the color lists given above, the list chromatic 



polynomial is 2. In contrast, for G = C 3 with color lists {ci} = (1,2,3), {c 2 } = (1,2), and 
{03} = (1) there is only one proper coloring, namely the color assignment (3, 2, 1) to vertices 
1, 2, 3, so the list coloring polynomial is equal to 1. 



XIII. CONCLUSIONS 

In this paper we have studied proper g-colorings of the vertices of a graph with a weighting 
factor w that either disfavors or favors a given color. In particular, we have analyzed a 
weighted chromatic polynomial Ph(G,q,w) associated with this problem, which generalizes 
the chromatic polynomial P(G,q). Since Ph(G,q,w) can be obtained as a special limit 
of the Potts model partition function in an external magnetic field, its study represents a 
fruitful confluence of statistical mechanics and mathematical graph theory. We have found a 
number of interesting properties of this weighted chromatic polynomial. Among others, we 
have shown how it encodes more information about the graph G, as shown by the fact that it 
is able to distinguish between certain graphs that yield the same chromatic polynomial. We 
have given formulas for Ph(G,q,w) for various families of graphs G, including line graphs, 
star graphs, complete graphs, and cyclic lattice strip graphs. For w G (0,1), Ph(G,q,w) 
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effectively interpolates between P{G, q) and P(G, q — 1). Using our results, we have discussed 
the zeros of Ph(G, q, w) in the q and w planes and their accumulation sets in the limit of 
infinitely many vertices of G. Finally, we have mentioned some observations, conjectures, and 
related weighted graph-coloring problems. There is ample motivation for further research on 
this very interesting subject. 
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TABLE II: Table of numbers np(L y ,d) and their sums, Np t L \ for the chromatic polynomial of cyclic 
strips of the lattice A (sq, tri, he) with h = 0. Blank entries are zero. See text for further discussion. 
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APPENDIX A: TABLES ON STRUCTURAL PROPERTIES 

For comparison with our new results for nph(L y , d) and Nph^ yl we list here corresponding 
tables for the following numbers for cyclic strips of the square (sq), triangular (tri), and 
honeycomb (he) lattices A: (i) np(L y , d) and their sums, Np t L x f° r the chromatic polynomial 
with h = 0, (ii) riz{L y , d) and their sums, Nz,l v ,x f° r the Potts model partition function with 



h = 0, and (hi 



with h + [J, [J. 



nzh{L y , d) and their sums, N Z h,L v ,\ for the Potts model partition function 
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TABLE III: Table of numbers nz(L y ,d) and their sums, Nz,g,x, for the Potts model partition function on 
cyclic strips of the lattice A (sq, tri, he) with h = 0. Blank entries are zero. See text for further discussion. 
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TABLE IV: Tabic of numbers nzh(L v , d) and their sums, Nzh,L y for the Potts model partition function on 
strips of the lattice A (sq, tri, he) with h ^ 0. Blank entries are zero. See text for further discussion. 
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APPENDIX B: Ph(G, q, w) FOR TREE GRAPHS G 



1. n = 5 Vertices 

There are three tree graphs with n = 5 vertices, as shown in Fig. Fig. [TJ (i) the line 
graph L 5 , (ii) the graph Y" 5 , and (hi) the star graph S5. The weighted chromatic polynomials 
for these are 



Ph(L 5 ,q,w) = (q-1) q 4 + (5w - 8)q 3 + 3(2w 2 - 9w + 8)q 2 



+ (w - 2){w 2 - lQw + 16)g - (w - l){w 2 - 12w + 16) 



q 5 - (9 - 5w)q 4 + 2(w - A)(3w - 4)g 3 - (-w 3 + 2Aw 2 - 75w + 56)g 2 



+ (-2w 3 + 31w 2 - 76w + 48)g - (1 - w){w 2 - 12w + 16) 



(Bl) 



Ph(Y 5 , q,w) = (q-l)(q + w- 2) q 3 + 2(2w - 3)q 2 + {2w 2 - 13w + 12) q - (w - l)(3w - 8) 



g 5 - (9 - 5w)q 4 + 2{w - 4)(3w - 4)g 3 - 2{-w 3 + 13w 2 - 38w + 28)q 2 



+ (-5w 3 + 37w 2 - 79w + A8)q - (w - l)(w - 2) (8 - 3w) 



(B2) 



Ph(S 5 ,q,w) = {q-1) q 4 + (5w - 8)g 3 + 3(2w 2 - 9w + 8)q 2 



+ (Aw 3 - 24w 2 + 51w - 32)g + (w - l)(w 3 - 7w 2 + 17w - 16) 



q 5 - (9 - 5w)q 4 + 2(w - A)(3w - 4)g 3 - 2(-2w 3 + 15w 2 - 39w + 28)q 2 



+ (w 4 - I2w 3 + A8w 2 - 8Aw + 48)g - {w - l)(w 3 - 7w 2 + 17 w - 16) (B3) 



2. Ph(G, q, w) for Tree Graphs with n = 6 Vertices 

There are six tree graphs with n = 6 vertices, as shown in Fig. [6j (i) the line graph Lq, 
(ii) the graph Ys, (iii) the graph with a branch in the middle of the line, denoted iso — Ye 
(iv) a graph with two branches, denoted H e , (v) a graph forming a cross, denoted CV 6 , and 
(vi) the star graph Sq. (Again, we order these in terms of graphs with increasing maximal 
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iso — Ya 




Ha 



Cra 



FIG. 6: Tree graphs with n = 6. 



vertex degree A; one has A = 2, 3, 3, 3, 4, 5 for graphs (i)-(vi), respectively. In chemical 
nomenclature, five of these graphs correspond to the carbon backbones of the following 
alkanes: (i) n-hexane, (ii) 2-methylpentane, (iii) 3-methylpentane, (iv) 2,3-dimethylbutane, 
and (v) 2,2-dimethylbutane. The graph Sq has no carbon-atom correspondence, since the 
central vertex has degree 5. We find 

Ph(L 6 , q, w) = (q- 1) \q 2 + 2{w - 2)q - 3w + 4] [g 3 + 2{2w - 3)q 2 



+ (2w 2 - 13w + 12)q - 2(w - l)(w - 4) 



q 6 - (11 - 6w)q 5 + 10(w 2 - 5w + 5)q 4 - 2(-2w 3 + 29w 2 - 82w + 60)q 3 



+ (-Uw 3 + 123w 2 - 264w + 160)g 2 - (-16w 3 + 113w 2 - 208w + 112)g 



+ 2(1 -w)(4-3w)(4- w) 



(B4) 



Ph(Y 6 ,q,w) = (q-l) q 5 + 2(3w - 5)g 4 + 2(5w 2 - 22w + 20)q 3 



+ (5w 3 - 50w 2 + 121w - 80)g 2 + (w 4 - I6w 3 + 8Aw 2 - U8w + 80)q 



(w - l)(w - 4)(w 2 - 7w + 



q 6 - (11 - 6w)q 5 + 10(w 2 - 5w + 5)q 4 - 5(-w 3 + 12w 2 - 33w + 24)q 3 
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+ (w 4 - 21w 3 + 134w 2 - 269w + 160)g 2 
- (2w 4 - 28w 3 + ISlw 2 - 216w + 112)g 

+ (w-l)(w-4)(w 2 -7w + 8) (B5) 

Ph(IsoY 6 ,q,w) = (q-l)(q-2 + w) q 4 + (5iu - 8)g 3 + (w - 4)(5w - Q)q 2 

+ (-14w 2 + A5w - 32)q + 2(w - l)(5w - 8) 

= g 6 - (11 - 6w)q 5 + 10(w 2 - 5w + 5)g 4 - 5(-u» 3 + 12w 2 - 33w + 24)g 3 

+ (-19w 3 + 133w 2 - 269w + 160)g 2 - (-24w 3 + 129w 2 - 216w + 112)g 

+ 2(w - l)(iu - 2)(8 - 5w) (B6) 

Ph(H 6 ,q,w) = (q- l)(g-2 + w) 2 g 3 + 2(2w - 3)g 2 + (w 2 - 12w + 12) q - 2(w - l)(w - 4) 

= g 6 - (11 - 6w)q 5 + 10(w 2 - 5w + 5)g 4 - 2(-3w 3 + 31w 2 - 83w + 60)g 3 

+ (w 4 - 26w 3 + UAw 2 - 27 'Aw + 160)g 2 - (w - 2)(3w 3 - 32w 2 + 8Aw - 56)g 

+ 2{w- l)(w -2) 2 (w -A) (B7) 

Ph(Cr 6 , q,w) = (q - l)(q - 2 + w) q 4 + (5w - 8)g 3 + (w - A)(5w - 6)q 2 
+ (2w 3 - 18w 2 + A7w - 32)q - (w - l)(3w 2 - 13w + 16)] 
= g 6 - (11 - 6w)q 5 + 10(w 2 - 5w + 5)g 4 - (-7w 3 + 6Aw 2 - 167w + 120)g 3 
+ (2w 4 - 32w 3 + 153w 2 - 278w + 160)g 2 - (5w 4 - A7w 3 + 160w 2 - 229w + 112)g 
+ (w - - 2)(3w 2 - 13w + 16) (B8) 

Ph(S 6 ,q,w) = (q-1) g 5 + 2(3w - 5)g 4 + 2(5w 2 - 22w + 20)g 3 
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+ 2(5w 3 - 30w 2 + 63w - 40)g 2 + (5w 4 - 40w 3 + 120w 2 - 164w + 80)q 



+ (w - l)(w 4 - 9w 3 + 31w 2 - A9w + 32) 



11 - 6w)q 5 + I0(w 2 - 5w + 5)g 4 - 10(-w 3 + 7w 2 - 17 w + I2)q 3 



+ 5(w 4 - 10w 3 + 36w 2 - 58w + 32)q 2 



-w 5 + 15w 4 - 80w 3 + 200w 2 - 245w + 112) q 



+ (1 - w)(w 4 - 9w 3 + 31w 2 - 49w + 32) 



(B9) 



Among other things, these calculations can be used to characterize further the way in 
which the weighted chromatic polynomial is able to distinguish between graphs that yield the 
same chromatic polynomial. As discussed in the text, all tree graphs with a given number n 
of vertices yield the same chromatic polynomial, P(G tre e,n, <?) = <?(<? — l)™^ 1 (and, indeed, also 



the same Tutte polynomial T(G tr ee n ,x,y) 



x 



n-l 



). Using our results above, we calculate 



the following differences in weighted chromatic polynomials, relative to Ph(Se,q,w), for 
defmiteness, from which all other differences can be obtained: 

Ph(S 6 , q, w) - P/i(L 6 , q, w) = w(w - l) 2 (q - 1) \(3q + w)(2q + w) - 20q - 8w + 17] (BIO) 

Ph(S 6 , q, w) - Ph(Y 6 , q, w) = w(w - lf(q - l)(5g 2 + Awq + w 2 - 16q - 7w + 13) (Bll) 

(B12) 



Ph(S 6 , q, w) - Ph(IsoY 6 , q, w) = w(w - l) 2 (q - 1) 
Ph(S 6 , q, w) - Ph(H 6 , q, w) 



5q 2 + 5wq + w 2 — I6q — 8w + 13 



w[w 



l) 2 (q-l)(2q-3 + w) 



3q 2 + 3wq + w 2 — 9q 



5w + 7 



(B13) 
(B14) 



Ph(S a , q, w) - Ph(Cr G , q, w) = w(w - lf(q - 1) 

We thus find that the weighted chromatic polynomials for all of the different n-vertex tree 
graphs of a given n are, in general, different from each other, although they coincide for 
w = 1 and w = 0, where they reduce to chromatic polynomials, and for 5 = 1, where they 
all vanish. 
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